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PART  I 
INTRODUCTION 


CHAPTER  1 
OVERVIEW 


Modern  society  is  more  and  more  revolutionized  by  technological  breakthroughs. 
Replacement  of  massive  manual  labor  by  machinery  has  been  completed  in  the  developed 
countries.  Great  progress  in  communication,  transportation,  energy  consumption,  health- 
related  industry,  and  so  on,  move  the  world  towards  a  more  environment-friendly,  long- 
lifespan  biosphere.  And  materials  research  and  development  (R&D)  has  tremendously 
facilitated  this  advancement.  For  instance,  steel  made  the  first  industrial  revolution  into 
a  reality;  copper  wire  enables  the  long  distance  transportation  of  electricity;  aluminium 
allows  airplanes  to  take  to  the  sky;  plastics  are  widely  used  in  packaging;  and  ceramics 
are  versatile  in  insulation.  But  the  semiconductor  presents  the  best  example. 

Semiconductors  have  been  making  resounding  impacts  upon  the  way  of  life  of 
mankind  since  Bardeen  and  Brattain  [I,  2]  discovered  the  point-contact  transistor,  and 
Shockley,  Sparks,  and  Teal  [3]  first  observed  the  junction  transistor  action.  This 
realization  of  the  solid-state  version  of  the  vacuum-tube  triode  began  the  era  of  circuit 
integration  and  a  computing  revolution.  The  rapid  progress  in  compound-semiconductor 
technology  introduces  a  new  class  of  ultrafast  devices — heterostructure  devices.  Without 
them  for  information  processing,  high-speed  computing  would  not  be  possible.  However, 
it  would  be  ignorant  without  mentioning  that  optoelectronic  devices,  which  are  essential 
in  solid  state  lasers,  telecommunication  and  display  technologies.  They  are  transforming 
the  field  of  telecommunications.  Nowadays,  semiconductors  research  is  indispensable  in 
the  success  and  development  of  the  microelectronic  and  the  optoelectronic  industries. 


3 
A  semiconductor  is  a  solid  with  an  electric  conductivity  between  that  of  a  metal 
and  an  insulator  [4j.  Discussion  in  this  overview  is  restricted  to  the  sub-group  of  IV, 
III-V,  and  II- VI  crystalline  compound  semiconductors.  Because  of  the  nature  of  the 
predominantly  covalent  bonding,  group  IV  elements  crystallize  in  diamond  lattice  (space 
group  07h);  most  III-V  compounds,  e.g.,  GaAs,  crystallize  in  zincblende  lattice  (space 
group  T$;  see  Fig  1),  many  of  the  II- VI  compounds  crystallize  in  both  the  zincblende 
and  the  wurtzite  lattices  (space  group  C$v)  [5],  as  tonicity  of  the  bonding  increases. 
Their  electric  conducting  ability  ranges  from  a  semi-metal  (e.g.,  a-Sn  and  HgTe)  to  a 
widegap  semiconductor  (e.g.,  CdS  and  GaN),  which  means  a  small  conduction-valence 
band  overlap  to  a  bandgap  as  iarge  as  6.2eVfor  AIN.  The  electronic  and  phonon  structures 
of  these  ciystalline  semiconductors  are  predominantly  determined  by  the  nature  of  their 
chemical  bonding  [6],  which,  in  turn,  determines  the  prospects  of  a  potential  candidate 
of  technical  and  scientific  interest 

The  unique  ability  for  semiconductors  to  be  modified  extensively  by  external  means 
is  the  very  reason  why  they  are  so  versatilely  employed  in  transistors,  switching  devices, 
voltage  regulators,  photocells,  photodetectors,  and  so  on.  Doping  is  the  original,  and 
still  is,  the  most  effective  way  to  modify  the  conducting  properties  of  a  semiconductor 
[7]  In  order  to  achieve  this,  shallow  impurity  states  are  desired  so  that,  under  thermal 
excitation,  free  carriers  can  be  generated  into  the  conduction  (electron  conduction)  and/or 
valence,  band  (hole  conduction).  Bipolar  and  p-n  junctions  are  made,  by  means  of  doping, 
for  circuit  switching,  signal  amplification,  and  current  rectifying.  Even  though  doping 
supplies  a  semiconductor  with  carriers,  it  restricts  its  mobility  by  presenting  scattering 
centers  in  the  system.  Modulation  doping  greatly  reduces  this  restriction  by  separating 
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5 
the  carriers  and  impurities  in  different  regions.  Over-doping  a  bulk  semiconductor  can 
convert  it  into  a  metal,  and  this  process  is  often  called  metallization.  On  the  other  hand, 
the  luminescent  properties  of  a  semiconductor  can  be  modified  by  doping  the  material 
with  deep-level  impurities,  since  the  quantum  efficiency  of  photon  generation  intricately 
depends  upon  the  interplay  of  radiative  and  non-radiative  channels. 

By  physically  and/or  chemically  controlling  the  interfacial  abruptness,  one  can  make 
semiconductor  heterostructures.  This  introduces  a  whole  new  world  of  opportunities  for 
new  concepts  for  devices  and  physics  [8-1 1  J.  Various  quantum  devices  are  made  into  re- 
ality with  progress  in  growth  technologies.  Quantum  wells,  wires,  dots,  and  superlattices 
allow  improved  electronic  properties  from  the  bulk  system,  enhanced  oscillator  strengths, 
and  density  of  states.  Wave-packet  dynamics  in  quantum  confined  systems  can  be  used  to' 
implement  Bloch  oscillation,  and  to  generate  Tera-Hertz  electromagnetic  waves.  Just  as 
carriers  (electrons  and  holes)  are  restricted  in  real  space,  confinement  may  also  occur  for 
phonons  and  photons.  In  heterostructures,  there  are  many  physical  parameters  which  can 
be  manipulated  for  the  advantage  of  devices.  Bandgap  engineering  [12],  which  means 
targeted  alteration  of  bandstructure,  can  be  done  by  alloying  in  the  barrier  and/or  active 
region,  introducing  strain  effect  due  to  lattice  mismatch,  or  creating  internal  electric  field 
due  to  charge  transfer.  Nanostructured  systems1  [13]  condense  many  developments  in  the 
area  of  heterostructures  and  play  an  irreplaceable  role  in  integrated  electro-machinery, 
electronics  and  optoelectronics.  Heterostructures  make  high-mobility,  high-efficiency, 
low-threshold,  'ow-dimension  possible  in  semiconductors. 

External  fields  can  manipulate  the  response  of  semiconductors  so  as  to  achieve  what 


1.  Low-dimensional  structures  that  exhibit  manifest  quantum  effects. 
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is  impossible  at  equilibrium.  High  electric  fields  can  modify  the  electronic  structures 
and  dynamics  of  semiconductors.  Two  dimensional  electron  gas  exhibits  novel  physical 
phenomena,  such  as  mesoscopic  transport,  quantum  Hall  effect,  and  enhanced  many- 
body  effects  [14].  Also,  under  high  electric  field,  hot  elections  are  generated,  due 
to  acceleration  by  the  field  and  bottlenecks  in  the  relaxation.  Some  systems  show 
the  Gunn  effect  [15]  in  bulk  semiconductors  and  avalanche  processes  in  reversed  p-n 
junctions  [16].  Thin-film  electric-field-related  displays  [17]  are  under  active  pursuit  in  the 
phosphor  community  for  high  electric  field  inside  a  semiconductor  thin  film  efficiently 
yields  visible  light,  under  certain  circumstances.  Magnetic  field  tends  to  localize  and 
polarize  electronic  states  with  Landau  levels  formed  in  a  semiconductor,  which  manifests 
interesting  magneto-optical  effects  [18]  and  quantum  Hall  effects  [19,  20].  Nowadays, 
powerful  lasers  are  not  only  useful  in  spectroscopy  [21]  but  also  capable  of  driving  a 
semiconductor  far  from  equilibrium.  Femtosecond  spectroscopy  provides  information  on 
the  dynamical  behavior  of  carriers.  New  physical  systems,  such  as  excitonic  plasma,  are 
artificially  created  to  test  theoretical  predictions,  such  as  Bose-Einstein  condensation  [22]. 
Though  some  important  methods  of  manipulation,  such  as  current  injection  and  sound 
waves,  are  neglected  in  the  present  discussion,  it  is  quite  amazing  that  in  so  many  ways 
external  sources  can  be  employed  to  detect,  modulate,  and  manipulate  a  semiconductor 
system. 

Above  is  an  incomplete  general  summary  of  semiconductor  physics  and  technology. 
Reliability  and  magnetic-effects-related  topics  are  not  covered.  In  the  next  chapter,  the 
author  will  focus  on  three  II- VI  semiconductors,  ZnS,  ZnSe  and  CdTe.  The  following 
diagram  illustrates  the  logical  connections  of  this  research. 
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CHAPTER  2 
IIB-VIA  COMPOUND  SEMICONDUCTORS 


The  present  thesis  deals  exclusively  with  zincblende,  widegap  chalcogenide  semicon- 
ductors. All  of  the  three  materials  that  are  considered  possess  a  direct  bandgap,  which 
is  instrumental  in  optical  and  optoelectronic  applications.  Narrow-bandgap  chalcogenide 
semiconductors  are  extensively  investigated  for  their  involvement  in  infrared  photoemis- 
sion,  detection  and  optical  modulation  [23].  The  need  for  high-speed,  wide-bandwidth 
transmission  of  information  pushes  for  R&D  into  the  widegap  semiconductors  [24],  both 
III-V  compounds  and  IIB-VIA  compounds  which  are  composed  of  Zn,  Cd,  Hg  (group 
1TB)  and  S,  Se,  Te  (group  VIA).  In  accord  with  the  progress  in  processing  technol- 
ogy for  these  compounds,  the  study  of  IIB-VIA  chalcogenides  inevitably  benefits  and 
consequently  Hg/CdTe  detectors  [25].,  ZnSe-based  blue/green  LEDs  [26],  and  other  op- 
toelectronic applications  [27]  are  available  commercially  today. 

ZnS 

As  one  of  the  most  important  phosphors  [30]  and  ultraviolet  (UV)  coating  materials 
[31],  great  efforts  in  academia  and  industry  have  been  dedicated  to  ZnS.  Theoretically, 
the  electronic,  transport,  and  optical  properties  of  ZnS  have  been  investigated  thoroughly 
(for  a  complete  list  of  references,  see  the  bibliography  in  Cohen  and  Chelikowsky 
[32]).  Nevertheless,  there  are  rejuvenating  activities  in  ZnS,  due  to  interest  in  the 
alternating-current  thin-film  electroluminescent  display  (AC  TFELD)  research  [33-36] 
and  the  coming  of  a  heterostructure  age. 
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Table  1:  Collection  of  Important  Parameters  of  ZnS,  ZnSe,  and  CdTe  [4,  28,  29]. 


Parameter 

ZnS 

ZnSe 

CdTe 

Lattice  constant,  a  (A) 

5.4093 

5.6687 

6.481 

Energy  gap  at  0  K,  Eg  (eV) 

3.34 

2.94 

1.59 

^  (10-4  eV/K) 

-4.6 

-8 

-3 

^f  (10"6eV/bar),  (at  R.T.) 

5.7 

0.7 

8.0 

Spin-orbit  splitting,  Ao  (eV) 

0.072 

0.403 

0.811 

Electron  effective  mass,  m*  (me) 

0.27 

0.17 

0.11 

Heavy  hole  effective  mass,  m\h 

0.58 

0.60 

0.40 

Light  hole  effective  mass,  m*lh 

0.230 

0.149 

0.13 

Valence  band  parameters: 

71 

6.12 

4.3 

5.3 

72 

0.53 

1.14 

1.7 

73 

1.00 

1.84 

2.0 

High  frequency  dielectric  constant,  e<x 

5.7 

6.2 

7.2 

Low  frequency  dielectric  constant,  cq 

8.9 

9.6 

10.4 

Electroluminescence  (EL)  is  a  cold  light-emitting  phenomenon  when  a  high  electric 
field  is  applied- onto  a  non-metallic  material.  With  the  advent  of  VLSI  technology, 
TFELD  becomes  a  viable  alternative,  among  flat  panel  displays  (FPD),  to  the  existing 
bulky,  high-voltage  cathode-ray  tube  (CRT)  technology.  There  presents  a  challenge, 
however,  to  the  R&D  in  phosphorescent  semiconductors,  which  is  how  to  improve  their 
performance  and  reduce  their  price.  In  technical  respect,  understanding  basic  physical 
processes,  such  as  carrier  dynamics,  radiative  relaxations  that  control  the  EL  phenomenon 
and  then  optimizing  TFEL  devices  are  the  major  concerns  of  display  R&D. 

To  understand  EL  in  a  Metal-Insulator-Semiconductor-Insulator-Metal  (MISIM) 
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structure  (see  Fig.  3),  free  carrier  dynamics  has  to  be  investigated  for  it  determines 
the  quantum  efficiency  of  photon  generation  by  impact  excitation  of  the  luminescent  im- 
purities (activators).  For  dilute  free  carrier  systems  that  lack  long  range  interactions,  a 
kinetic  description  of  the  carrier  distribution  functions  has  been  found  to  be  overwhelm- 
ingly successful  [37-42].  However,  Boltzmann  kinetic  equation  (BKE)1  is  notoriously 
difficult  to  solve  analytically.  Nowadays,  solutions  of  BKE  require  numerical  compu- 
tations and  Monte  Carlo  simulations  [43-45].  These  require  immense  computing  time 
and  restrict  their  applications  to  large-scale  dynamics  modeling.  Full  Brillouin  zone  elec- 
tronic structures  that  are  required  in  sophisticated  computations  increase  the  demand  for 
computing  power  considerably. 

Inhomogeneity  in  the  MISIM  EL  system  introduces  spatial  degrees  of  freedom  into 
BKE  and  demands  self-consistent  treatment  with  respect  to  the  non-uniform  local  electric 
field.  To  meet  this  numerical  challenge,  numerical  integration  scheme,  i.e.,  adaptive 
stepsize  Runge-Kutta  method  [46,  47],  may  be  adopted  to  solve  BKE  with  the  electronic 
structure  of  the  system,  derived  from  the  fast  empirical  30  x  30  k-p  method  (see  Appendix 
A),  included.  Part  II  of  the  present  thesis  covers  preparatory  research  for  this  future 
project.  Bandstructure  calculations  for  ZnS  are  implemented  and  employed  for  modeling 
the  X-ray  photoemission  spectroscopy  spectrum  and  optical  reflectivity  measurement. 
Simultaneously,  the  calculated  bandstructure  is  compared  to  well-established  experimental 
and  theoretical  results  for  accuracy  check.  The  results  are  in  excellent  agreement  with 
the  well-established  data  and  fit  the  measurements  satisfactorily. 


1 .  also  known  as  Boltzmann  transport  equation  (BTE). 


12 
ZnSe 


Due  to  its  similarity  to  ZnS,  research  in  ZnSe  has  always  been  done  in  conjunction 
with  the  research  in  ZnS  [28].  The  major  differences  between  these  two  materials  are 
twofold:  (1)  ZnS  has  an  energy  gap  of  3.84  eV  and  ZnSe  possesses  one  of  2.94  eV;  (2) 
spin-orbital  coupling  in  ZnSe  is  much  more  important  than  in  ZnS,  thanks  to  the  lar«e 
atomic  weight  of  selenium.  However,  their  technological  importance  is  more  related  to  the 
availability  of  high  quality  materials  that  industry  can  produce.  In  applications  concerning 
phosphors,  the  materials  are  in  polycrystalline  structures  while  modern  developments  are 
more  and  more  inclined  towards  single  crystalline  structures.  In  this  aspect,  ZnSe  is 
chosen  more  than  ZnS  in  quantum  structure  R&D.  It  is  not  within  the  scope  of  this 
dissertation  to  provide  a  review  of  bulk  ZnSe.  For  a  list  of  references  on  the  physical 
properties  of  bulk  ZnSe,  see  the  bibliography  in  Cohen  and  Chelikowsky  [32]. 

The  current  research  interests  in  ZnSe  concentrate  upon  the  area  of  blue/green  LEDs 
and  lasers  [26].  These  quantum-well  devices  with  ZnSe,  as  the  confined  and  active 
region,  sandwiched  between  barrier  layers,  provide  light  as  a  result  of  carrier  injection 
through  the  electrodes  and  radiative  recombination  in  the  confined  region.  One  of  the 
technical  challenges,  once  severely  impeding  R&D  in  these  devices,  is  low-resistivity 
p-type  shallow-impurity  doping  in  ZnSe,  which  is  required  for  making  ohmic  contacts 
in  injection-type  LEDs  and  lasers.  It  was  solved  by  the  ingenious  method  of  Park  et 
al  [48].  Even  though  available  on  the  market  today,  there  are  still  plenty  of  interesting 
physics  to  be  done  in  these  quantum-structured  devices.  Among  them  is  the  study  of  the 
physical  properties  of  intermediately  to  heavily  doped  bulk  systems.  Optical  properties 
of  these  systems  affect  the  performance  of  the  doped  semiconductors  as  the  transparent 
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window  for  light  output.  Transport  properties  are  the  major  issue  in  the  making  of  ohmic 
contacts  and  quantum  efficiency  of  these  injection-type  devices. 

Part  IV  of  the  present  thesis  treats  the  «-type  shallow-impurity  doping  problem  in 
ZnSe.  The  major  difference  in  treating  p-type  doping  is  that  the  transport  is  hole-type 
and  degeneracy  in  the  bands  complicates  the  problem.  The  n-type  doping  problem 
has  been  treated  by  Harry  Ruda  [49].  A  shortcoming  of  his  analysis  is  it  does  not 
consider  the  change  in  the  chemical  potential  when  the  system  undergoes  intermediate  to 
heavy  doping.  In  reality,  there  are  different  kinds  of  extrinsic  impurities,  in  addition  to 
native  point  defects  and  their  complexes  in  the  rc-type  doped  ZnSe.  The  present  research 
considers  the  most  simplified  case,  which  is  a  system  with  only  independent  hydrogen- 
like dopants  and  acceptors  under  the  assumption  of  complete  compensation,  but  takes 
into  account  self-consistent  static  screening  for  impurity  ionization.  The  model  allows 
the  chemical  potential  to  change  with  doping  and  thus  modify  the  carrier  concentration 
as  well  as  the  mobility  correspondingly.  Qualitative  agreements  for  chemical  potential, 
carrier  concentration  are  found  between  the  related  experiment  and  the  present  research, 
but  lack  of  low  temperature  measurements  on  the  mobility  impedes  comparison  in  this 
aspect. 

CdTe 

Crystalline  cadmium  telluride  [50]  has  been  employed  in  nuclear  radiation  detection, 
electro-optic  modulation,  and  as  an  infrared  optical  material.  Alloys  of  CdxHgi_xTe  have 
become  important  as  infrared  detectors,  and  CdTe  is  a  key  substrate  for  the  fabrication 
of  devices  based  on  CdxHg[_xTe  alloys  and  superlattices.  Potential  applications  in  EL 
and  microwave  devices,  high  power  laser  windows,  etc.,  have  stimulated  considerable 
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interests  of  basic  research  in  CdTe.  Optical  and  photoemission  studies  were  undertaken  in 
aid  of  the  construction  of  CdTe  bandstructure  (  see  reference  50,  page  84).  The  electronic 
structure  of  CdTe  was  studied  by  applications  of  the  dielectric  bonding  theory  of  Philips 
and  Van  Vechten  [6]  and  various  other  more  detailed  approaches  (  see  reference  50, 
page  81),  which  provide  good  interpretations  to  experimental  data.  Transport  properties 
of  CdTe  provide  further  information  such  as  the  effective  masses  of  charged  carriers, 
interactions  of  charged  carriers  and  phonons  with  impurities  and  elementary  excitations 
in  the  material.  Moreover,  defects  in  CdTe  are  extensively  investigated  because  their 
structures  influence  performances  of  CdTe  in  devices.  For  the  latest  experimental  and 
theoretical  data,  Landolt-Bornstein  New  Series  [29],  pages  225-230,  may  be  consulted. 

Semiconductor  microstructures  and  nanostructures  are  important  in  integrated  systems 
and  broaden  the  possibilities  of  improving  optoelectronic  properties  of  semiconductors 
by  exhibiting  behavioral  transitions  from  bulk  materials  to  molecules.  These  quantum 
size  effects  in  CdTe  micro-  and  nanostructures  have  been  explored  extensively  in  the  last 
twenty  years  [13,  28].  Quantum  wells,  wires,  dots,  and  superlattices  are  configurations 
of  those  well-studied  quantum  structures.  Among  these  structures,  quantum  dots  (QDs) 
show  carrier  confinement  effects  in  all  three  dimensions,  which  receive  a  great  deal  of 
attention  for  their  significant  optical  nonlinearity.  One  way  to  make  QDs  is  chemically 
embedding  isolated  semiconductor  crystallites  in  a  transparent  glass  matrix  [51],  and 
these  systems  are  widely  investigated.  CdS  and  CdTe  crystallites  are  the  main  focus  of 
researches  because  of  their  utility  in  integrated  optical  systems.  Systematic  researches 
of  the  quantum  confinement  effects  (QCEs)  on  various  properties  of  nanostructures  have 
been  carried  out  to  elucidate  fundamental  physics  in  these  systems  [42,  52-57],  not 
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necessarily  only  on  CdTe. 

Despite  impressive  efforts,  there  are  issues  in  CdTe  crystallites  that  remain.  One  of 
them  is  the  QCEs  on  the  absorption  edge.  The  saturation  of  QCEs  at  extreme  confinement 
seen  by  Potter  etal  [58]  was  shown  to  disagree  with  current  theoretical  models.  To  better 
model  QCEs  in  widegap  QDs,  it  is  essential  to  realize  that  although  Coulomb  interaction  is 
important  in  modifying  the  oscillator  strengths  of  optical  transitions,  bandstructure  effect 
is  predominant  in  affecting  the  optical  transition  states.  Part  IE  of  the  present  thesis 
deals  with  the  problem  along  this  line  with  an  approximate  approach — the  truncated 
crystal  approximation,  with  the  bulk  bandstructure  given  by  the  empirical  30  x  30  k  ■  p 
method.  Excellent  agreement  is  found  between  the  experiment  of  Potter  et  al.  and  the 
present  theoretical  work. 

As  a  note,  CdTe  is  classified  in  this  research  as  a  widegap  semiconductor,  despite 
the  fact  that  it  has  a  bandgap  of  1.59  eV.  This  is  simply  because  its  bandstructure  (Fig. 
13)  greatly  resembles  those  of  the  other  widegap  semiconductors,  such  as  ZnS  (Fig.  4). 
As  far  as  the  study  of  its  optical  properties  is  concerned,  no  differences  are  discerned 
compared  with  the  other  widegap  systems  in  reality. 


PART  II 
OPTICAL  PROPERTIES  AND  XPS  SPECTRUM  OF  ZnS 


CHAPTER  3 
INTRODUCTION 


Zinc  sulfide  has  been  one  of  the  most  important  phosphors,  used  as  visual  displays 
in  scientific  instruments,  entertainment,  and  radars.  It  is  also  an  indispensable  coating 
material  in  applications  of  ultraviolet  optics.  Though  tremendous  progress  has  been 
made  in  the  field  of  information  exchange  since  the  invention  of  optical  fiber,  high 
demand  on  speed  and  bandwidth  in  communications  has  always  been  in  place.  Faster, 
higher  frequency  and  more  efficient  optoelectronic  materials  are  under  active  pursuit.  In 
addition,  nonlinear  optical  properties  are  of  great  interest  in  telecommunications,  which 
is  instrumental  in  signal  amplification  and  modulation. 

Despite  the  fact  that  it  has  been  known  for  decades  as  a  phosphorescent  material, 
ZnS  has  not  finished  showing  improvements  in  efficiency,  spectrum,  and  refreshing 
rate.  For  the  last  three  decades,  great  interest  has  been  shown  in  development  of  the 
material  as  a  promising  active  blue/green  phosphor  [28].  Zincblende  ZnS  is  a  direct 
bandgap  semiconductor  with  a  bandgap  of  3.7  eV  at  room  temperature.  This  absorption- 
free  window  for  visible  light  and  excellent  properties  as  a  phosphor  make  ZnS  one  of 
the  top  candidates  for  R&D  in  solid  state  blue/green  light  emitting  diodes  (LEDs)  and 
lasers.  Study  of  optical  properties,  both  linear  and  nonlinear,  aids  understandings  towards 
development  and  optimization  of  microelectronic  and  optoelectronic  devices. 

In  the  present  part  of  this  research,  numerical  computation  of  the  linear  optical 
properties  from  microscopic  theory  is  needed  for  checking  bandstructure,  matrix  elements 
calculation,  and  Gilat-Raubenheimer  (GR)  scheme  code.    In  this  calculation,  (1)  the 
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empirical  30  x  30  k  p  full  zone  bandstructure  calculation  (see  Appendix  A)  for  zincbiende 
ZnS  (see  Fig.  4  or.  page  26)  has  been  implemented;  (2)  density  of  states  (DOS)  (see 
Fig.  5  on  page  28),  joint  density  of  states  (JDOS),  and  dielectric  functions  (see  Fig. 
8  on  page  32)  are  computed,  (3)  valence  band  DOS  of  ZnS  data  was  compared  with 
XPS  experiments  [59]  by  Dr.  Holloway  of  Materials  Science  and  Engineering  at  UF;  (4) 
calculated  dielectric  functions  are  examined  against  both  well-established  theoretical  and 
experimental  data  for  ZnS;  (5)  reflectivity  and  absorption,  contributed  by  linear  interband 
transitions,  are  computed  and  compared  to  experimental  data  by  Dr.  Holloway,  which 
show  good  agreement  (see  Fig.   9  on  page  34). 


CHAPTER  4 
THEORIES 

X-ray  Photoemission  Spectroscopy  (XPS)  [60] 

Photoelectron  spectroscopies  are  among  the  most  powerful  experimental  tools  for 
probing  the  electronic  structure  of  solids.  Photoemission  spectra  can  extract  information 
on  the  absolute  position  of  an  energy  band.  It  is  even  possible  to  determine  directly  the 
electronic  structure.  Also  photoemission  is  sensitive  to  surface  conditions,  either  intrinsic 
or  extrinsic  in  nature.  Careful  tuning  of  the  incident  radiation  enables  emphasizing 
either  bulk  or  surface  spectral  features.  Two  major  technical  developments  have  aided 
photoemission  measurements  immensely:  ultrahigh  vacuum  and  synchrotron  radiation. 
Photoemission  techniques  can  be  classified  according  to  how  the  energy  and  momentum 
of  the  incident  photon  and  ejected  photoelectron  are  controlled. 

The  three-step  model  is  a  most  often  used  physical  model  of  photoemission.  The  first 

step  is  a  photon  absorption  by  a  bulk  excitation.  Then  the  excited  electron  propagates 

to  the  surface,  and  the  third  step  is  its  escape  through  the  crystalline  surface.  Owing  to 

the  bulk  nature  of  the  first  step,  the  initial  and  the  final  wavevectors  are  identical  so  that 

a  direct  transition  occurs.  The  probability  for  an  excitation  to  occur  between  two  bands 

is  dictated  by  dipole  interactions,  and  energy  is  conserved  in  the  process.  The  number 

of  photoelectrons  ejected  with  a  specific  energy  for  a  given  incident  photon  energy  can 

be  formulated  as 

i(Ekin,hu)  =  Cl  x  J2  Piffysfafy  -  Ei(ic\  -  he 

xs(Ef{k)  -Ekin  +  Ev) 


i,i,f  '  (1) 


19 


20 
where  fiu  is  the  photon  energy,  Ekm  is  the  kinetic  energy  of  the  photoelectron,  Ev  is 
the  vacuum  level,  E, (k)  is  the  energy  of  the  initial  slate,  and  E}  \k\  is  the  energy  of 
the  final  state.  Pifikj  is  the  dipoie  matrix  element  between  initial  state  and  final  state. 
The  first  8  function  in  Eq.  (1)  comes  from  energy  conservation  in  the  first  step;  and  the 
second  8  function  imposes  energy  conservation  in  the  third  step. 

Under  certain  assumptions,  Eq.  (1)  reduces  to  an  extremely  useful  form.  For  instance, 
consider  the  assumptions  of  (1)  constant  matrix  elements  independent  of  k  and  of  (2)  flat 
valence  (or  conduction)  bands  or  non-conservation  of  k.  The  former  assumption  is  usually 
valid  though  the  latter  assumption  is  usually  not.  Under  these  additional  approximations, 
it  is  possible  to  obtain  [32] 

n(Ektn,hu)  =  c2  x  Ni(Ekin  -  Ev  -  hu) 

(2) 

*Nf(Ekin-Ev) 

where  c2  is  a  constant,  N,  is  the  DOS  for  the  initial  states,  and  Nf  is  the  DOS  for  the 
final  states. 

For  XPS,  the  mode  of  operation  is  to  measure,  as  a  function  of  energy,  the  number 
of  emitted  photoelectrons  "by  fixing  the  incident  photon  energy.  At  high  photon  energies, 
photoemission  spectra  are  fairly  easy  to  interpret.  From  an  operational  point  of  view, 
they  appear  to  replicate  the  valence  band  DOS.  In  a  more  formal  sense,  if  the  photon 
wavevector  is  not  insignificant  compared  to  a  reciprocal  lattice  vector,  which  may  be  the 
case  at  large  photon  energies,  then  ^-conservation  still  holds,  but  the  transitions  involved 
are  not  direct.  Therefore,  it  effectively  averages  over  the  final  states  in  Eq.  (1).  Moreover, 
at  large  photon  energies,  the  final  states  tend  to  be  plane-wave  like  and  most  likely  do  not 
contribute  to  a  structural  spectrum.  In  this  sense,  Eq.  (2)  becomes  valid.  Furthermore, 


21 
the  final  DOS  may  be  omitted  from  consideration: 

n(Ekin,  M  =  c  x  A^i(^iif,  -  Ev  -  ftw) .  (3) 

Therefore,   By  a  judicious  choice  of  hu>  and  via  an  analysis  of  Nu  a  direct  replica  of  the 
valence  band  density  of  states  is  almost  obtained    [32]. 

XPS  measures  energy  distribution  curves  (EDCs)  and  many  features  of  EDCs  are  well 
understood.  When  ftw  <  15eV,  EDCs  are  highly  structured,  which  arises  from  JDOS 
effects.  Other  "understood"  features  include  matrix  element  effects  and  penetration  depth 
changes.  The  matrix  element  effects  are  characteristic  of  the  electronic  structure  of  a 
semiconductor.  Above  the  fundamental  energy  gap.  dipole  transitions  occur  between 
upper  valence  bands  and  lower  conduction  bands.  The  initial  and  final  states  near  the 
band  gap  are  usually  well  defined  in  terms  of  orbital  character,  so  spectral  features  could 
be  enhanced  or  suppressed  depending  upon  the  magnitude  of  the  dipole  matrix  element 
between  valence  and  conduction  states.  At  higher  energies,  sum  rules  for  the  dipole 
matrix  elements  demand  that  the  magnitude  of  the  matrix  element  between  initial  and 
final  states  decrease.  As  a  consequence,  with  higher  energy  photons,  the  spectral  features 
begin  to  resemble  a  simple  JDOS.  The  critical  point  behaviors  in  the  spectral  features  are 
discussed  thoroughly  by  Cohen  and  Chelikowsky  [32]  and  Bassani  and  Parravicini  [61]. 

Linear  Optical  Properties — Interband  Transition  Theory  [61] 

The  optical  behavior  of  a  material  is  completely  determined  if  its  optical  constants 
are  known.  The  index  of  refraction,  n,  and  the  extinction  coefficient,  jfc,  may  thus  be 
regarded  as  a  summary  of  the  experimental  results  of  optical  measurements  made  on  the 
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material.  In  general,  they  are  related  to  the  dielectric  constants  by 

N2  =  e\  +  ie2 ; 


(4) 


N  =  n  +  iib , 

where  Ar  is  the  complex  index  of  refraction,  e\  and  e2  are  the  real  and  imaginary  parts 

of  the  complex  dielectric  constant  c.  The  normal-incidence  reflectivity  R  is  given  by 


R 


N-i" 


N  +  l 


n-lV  +  k2 

—o. -•  (5) 


In  this  part  of  the  thesis,  the  complex  dielectric  constant  t  is  computed  microscopically 
as  a  function  of  the  energy  of  incident  photon  hu  by  random-phase  approximation  (RPA) 
as  follows 

e(o,w\  =  1  +  4ttx(o,w) 

Im[x(6,u)}  =  ^Y,  J<*F(k)  6[u-uc,v(k)} 


cv  BZ 


Re 


cv  gz       UJ  -  uc%v  [k\ 


(6) 


F(k)  =  -^\(ck\e-P\uk)\2  . 
where  x  is  the  linear  response  function  of  the  system  to  the  transverse  electric  field  of 
the  incident  photon  with  polarization  e,  uCiV(k)  is  the  transition  frequency  between 
conduction  and  valence  bands,  and  jV  b  (^)2.  The  bands  are  computed  by  the 
empirical  30x30  k-p  method  (see  Appendix  A  and  references  within),  and  the  integration 
over  the  full  Brillouin  zone  (BZ)  is  implemented  via  the  Gilat-Raubenheimer  scheme  as 
outlined  in  Appendix  B  (see  references  within). 

Experimentally,  the  optical  properties  of  solids  may  be  determined  by  Kramers- 
Kronig  analysis  of  normal-incidence  reflectivity  data  or  indirectly  by  ellipsometric  mea- 
surements with  polarized  light.  It  is  convenient  to  distinguish  three  spectral  regions  in 
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the  description  of  the  optical  properties  of  widegap  semiconductor  compounds  The  first 
region,  extending  to  about  8  to  JO  eV,  is  characterized  by  a  sharp  structure  associated 
with  valence-  to  conduction-band  transitions.  The  second  region,  which  extends  to  about 
16  eV,  is  marked  by  a  rapid  decrease  of  the  reflectivity  that  is  reminiscent  of  the  behavior 
of  certain  metals  in  the  ultraviolet  Sharp  maxima  in  the  function  —Irn (e~~l),  which 
describes  the  energy  loss  of  fast  electrons  traversing  the  material,  have  been  frequently 
associated  with  the  existence  of  plasma  oscillation  in  this  region.  A  maximum  of  this 
type  is  seen  in  this  "metallic"  region.  In  the  third  region,  the  reflectivity  again  rises, 
indicating  the  onset  of  other  optical  absorption  processes.  This  structure  is  associated 
with  transitions  between  filled  d  bands  lying  below  the  filled  valence  bands  and  empty 
conduction  bands. 


CHAPTER  5 
RESULTS  AND  DISCUSSION 

The  electronic  structure  of  zincblende  ZnS  is  obtained  by  application  of  the  empirical 
30  x  30  k  ■  p  method  and  the  result  is  shown  in  Fig.  4.  Parameters  used  in  the  method 
are  adjusted  to  give  the  best  overall  fit  to  the  known  experimental  data  [29].  Tables  3 
and  2  list  the  k  ■  p  input  parameters  and  comparison  of  the  theoretical  result  to  well- 
established  experimental  data.  Refer  to  Appendix  A  for  an  explanation  of  the  meaning 
of  all  symbols.  V6~  and  E(rl25,)  are  set  to  zero,  so  the  top  of  the  valence  band  is  taken 
as  a  reference  point  for  the  band  energies. 

Table  2:  k  ■  p  Results  for  ZnS  Compared  with  Experimental  Values  [29]. 


Eg 

1  15  ~  l  15 

A0 

A3  -  Ai 

^3  ""  L3 

xi  -  xi 

expt. 

3.85 

8.35 

0.072 

5.7 

9.8 

6.6 

theory 

3.84 

8.76 

0.069 

5.95 

9.4 

6.88 

The  DOS  of  the  zincblende  ZnS  is  computed  as  well.  The  DOS  is  given  by 


N(E)  =  V^2 


d3k 


v2x)* 

71   BZ 


8  E 


*(*)) 


(7) 


where  n  is  the  band  index  which  takes  into  account  of  the  spin  degree  of  freedom  as 
well.  GR  scheme  is  used  when  calculating  the  integral  over  the  BZ.  Then  the  above 
expression  is  reduced  to 

V 


V       /        n     itetr=l     i      ■ 


(8) 
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Table  3:  30  x  30  k  ■  p  Method  Input  Parameters  for  ZnS,  ZnSe,  and  CdTe. 


Parameter 

ZnS 

ZnSe 

CdTe 

E(T[) 

-0.842 

-0.905 

-0.728 

E(TV 

0.943 

0.521 

0.6263 

E(T'2,) 

0.032 

0.0467 

0.078 

E(Tl) 

1.161 

1.028 

0.693 

E(TU.) 

0.755 

0.578 

0.6093 

E(T1S) 

0.433 

0.338 

0.416 

W^) 

1.1 

0.883 

0.8764 

p 

1.371 

1.056 

1.072 

Q 

1.16 

1.244 

1.035 

R 

0.847 

0.692 

0.769 

P" 

-0.014 

0.109 

0.157 

P' 

0.093 

0.232 

0.236 

Q' 

-0.691 

-0.508 

-0.526 

R' 

1.23 

1.355 

0.947 

pin 

1.31 

1.496 

1.021 

T 

1.379 

1.324 

1.034 

r 

-0.361 

-0.255 

-0.348 

x5 

0.006 

0.0179 

0.0542 

a15 

0.0013 

0.0245 

0.0409 

A- 

0.0 

-0.0141 

-0.0351 

Yf 

0.269 

0.3024 

0.1002 

v; 

•0.594 

-0.364 

-0.1503 

Vs~ 

0.76 

0.474 

0.260 

vr 

0.19 

0.3915 

0.12 

K 

-0.196 

-0.3043 

-0.0303 
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where  the  notation  follows  Appendix  B  and  the  integral  is  carried  out  inside  the  i-th 
smaller  tetrahedron.  Figure  5  is  computed  using  123  smaller  tetrahedra  within  each  of 
these  three  main  tetrahedra,  FLKW.  TLWU,  and  ITJWX.  The  computation  on  a  DEC 
alpha  3000/700  machine  takes  about  6  minutes.  Actually  the  number  of  smaller  tetrahedra 
above  may  be  reduced  to  63  without  too  much  sacrifice  of  numerical  accuracy.  The  total 
DOS  for  each  band  is  obtained  to  check  whether  it  agrees  with  unity,  the  expected  number 
of  electrons  per  unit  cell  that  each  band  can  take. 

Shown  in  Fig.  6  is  the  XPS  valence  band  data  after  removal  of  the  background 
and  shifting  1.1  eV  accordingly  to  fit  the  theoretical  data.  Two  subpeaks  (I,  and  I2)  are 
resolved  within  Peak  I  and  one  shoulder  (Si)  is  distinguished  in  the  spectrum.  In  order 
to  account  for  the  instrumental  broadening,  the  theoretical  DOS  data  is  convoluted  with 
a  Gaussian  distribution  with  a  full  width  at  half  maximum  (FWHM)  of  1.25  eV,  which  is 
the  FWHM  of  the  Zn  3d  core  level  in  the  XPS  spectrum.  The  result  is  shown  in  Fig.  7. 
It  is  easy  to  see  that:  (1)  the  relative  intensities  of  Peaks  I  and  II  fit  better  with  the  XPS 
data;  (2)  the  tail  at  the  high  energy  side  of  the  convoluted  curve  extends  to  about  1  eV 
above  the  original  top  of  the  valence  band.  Therefore,  the  theoretical  DOS  consistently 
explains  the  XPS  spectrum. 

The  real  and  imaginary  dielectric  constants  are  calculated  based  upon  the  bandstruc- 
ture  derived  by  the  empirical  30  x  30  k-p  method  before  other  optical  properties  can  be 
obtained.  The  calculation  is  carried  out  with  GR  scheme  and  an  approximation  of  constant 
optical  transition  probability  within  each  smaller  tetrahedron  is  used.  To  ascertain  the 
accuracy  of  the  calculated  dielectric  constants,  comparison  of  the  present  theoretical  calcu- 
lation with  well-established  experimental  data  (see,  for  example,  Landolt-Bbmstein  New 
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Series  [29],  pages  61-126)  is  made.  The  static  dielectric  constant,  e(0),  for  zincbiende 
ZnS  is  measured  to  be  from  8.0  to  8.9  at  room  temperature  and  this  research  finds  8.34  at 
3  eV;  the  optical  dielectric  constant,  e(oo),  is  measured  to  be  5  7  to  5.7  at  room  temper- 
ature and  this  research  obtains  5.20  at  10.5  eV.  (Normally,  for  widegap  semiconductors, 
the  relative  temperature  coefficient,  defined  as  —*§#-,  is  about  3  x  lO"4.  Therefore,  an 
increase  of  0.25  in  dielectric  constants  at  low  temperature  from  room-temperature  value 
is  expected.)  However,  this  comparison  has  to  be  taken  cautiously  because  the  optical  di- 
electric constant  does  not  originate  from  interband  transition  but  rather  from  high-energy 
collective  excitations  (plasmons)  of  the  material.  Therefore,  discrepancies  of  the  optical 
constants  from  the  linear  interband  transition  contributions  are  expected. 

As  a  note,  this  calculation  can  be  easily  extended  to  other  linear  response  functions. 

In  general,  a  linear  response  function  may  be  expressed  as 

1   ^         F[k) 
G(u)  =  lim  ~  Y }-l (9) 

if   «•-■«(*)  -*e 
The  complex  dielectric  function  is  a  special  case.  The  extension  is  straight  forward  given 

the  knowledge  of  F(k). 

-The  comparison  of  the  theoretical  reflectivity  is  done  with  the  experiments  carried 
out  in  Dr.  Holloway's  group  in  the  Department  of  Materials  Science  and  Engineering  at 
the  University  of  Florida.  This  calculation  assumes  zero  temperature  but  the  experiments 
were  done  at  room  temperature.  It  is  expected  that  details  that  are  prominent  at  low 
temperature,  such  as  critical  points,  tend  to  smear  out  at  room  temperature.  In  addition, 
energy  shifts,  as  much  as  200  rneV,  can  occur  due  to  temperature  effects.  In  Fig.  9, 
the  low  energy  oscillation  peak  in  experiment,  below  3.7  eV,  is  due  to  Fabry-Perot 
effect  in  the  thin-film  configuration  and  high  energy  data  ends  at  9.5  eV  due  to  shortage 
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of  tunable  light  source  and  lack  of  interest.  The  approximate  5%  difference  between 
experiments  and  theory  is  not  unusual.  A  possible  explanation  is  that  the  surface  of 
a  thin  film  sample  is  contaminated.  This  reduces  the  effectiveness  of  the  thin  film  to 
reflect  incident  light,  and  the  well-established  result  (see  Landolt-Bornstein  New  Series 
[29],  page  380)  tends  to  support  this  possibility.  Also,  the  approximation  of  constant 
optical  transition  probability  within  each  of  those  smaller  tetrahedra  introduces  further 
uncertainty  though  this  uncertainty  diminishes  when  the  number  of  partitions  for  the 
irreducible  BZ  is  increased.  A  check,  however,  showed  that  the  theoretical  reflectivity 
does  not  sensitively  change  with  the  number  of  partitions  of  the  irreducible  BZ. 
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CHAPTER  6 
CONCLUSIONS 

In  this  part  of  the  research,  the  X-ray  photoemission  spectrum  and  the  optical 
properties  of  zincblende  ZnS  are  theoretically  investigated.  The  empirical  30  x  30  k  ■  p 
method  has  been  employed  to  obtain  the  electronic  structure  in  full  BZ.  A  GR  scheme  was 
adopted  for  numerical  integration  over  the  BZ  to  find  both  DOS  and  dielectric  constants 
of  the  material.  Comparisons  with  well-established  experimental  data  were  made  to 
establish  the  accuracy  of  the  theoretical  method  and  also  allow  the  modeling  of  new 
experimental  data  by  Dr.  Holloway.  As  a  result,  the  following  conclusions  are  drawn: 

•  Bandstructure  calculation  is  realistically  accurate  in  the  full  Brillouin  zone; 
Up  to  10  compound  semiconductors  of  zincblende  structure  have  be  computed; 

•  Gilat-Raubenheimer  scheme  is  successfully  utilized  in  the  full  BZ  integration; 

•  Both  experimental  XPS  and  optical  properties  of  zincblende  ZnS  are  satisfactorily 
modeled. 
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PART  III 
QUANTUM  CONFINEMENT  EFFECTS  !N  CdTe 


CHAPTER  7 
INTRODUCTION 


Quantum  confinement  effects  (QCEs)  have  been  shown  to  cause  oscillatory  behaviors 
in  electronic,  thermodynamic  properties  for  metallic  [62,  63],  semimetallic,  and  semicon- 
ducting films  [64-66].  In  the  last  two  decades,  QCEs  have  been  widely  used  to  modify 
the  optoelectronic  properties  of  quantum  semiconductor  structures.  Electronic  properties 
which  depend  upon  QCEs  in  semiconductors  are  electronic  bandgap,  density  of  states, 
and  enhanced  oscillator  strengths  of  optical  transitions.  Different  configurations,  i.e., 
dot,  wire,  film,  and  superlattice,  are  adopted  to  implement  the  corresponding  design. 
The  improved  electronic  structures  allow  better  performances  of  the  related  quantum 
semiconductor  devices.  For  instance,  low  threshold  lasers  have  been  implemented  for 
GaAs/AlGaAs  quantum  well  (QW)  structures  [8]. 

Theoretical  efforts  in  the  explanation  of  QCEs  in  semiconductors  are  not  lacking[65, 
67,  68].  On  one  hand,  pseudopotential  method  and  fight-binding  method  take  into 
account  of  the  full  potential,  which  demands  large-scale  numerical  calculations.  On  the 
other  hand,  among  approximate  methods,  the  envelop  function  approximation  (ERA)  [69, 
70]  is  widely  used,  owing  to  its  simplified  factorization  of  the  electronic  wavefunction. 
However,  the  ERA  method  fails  in  strong  confinement  regime.  In  contrast,  the  truncated 
crystal  approximation  (TCA)  [71]  treats  the  enhanced  QCEs  in  thin  silicon  quantum 
films  successfully  [68]. 

The  TCA  method  predicts  the  electronic  properties  of  the  quantum  semiconductor 
structures  from  those  of  the  bulk  materials,  and  vice  versa.     Quantum  confinement 
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effects  in  semiconductors  are  treated  under  the  TCA  method  by  applying  semiclassica! 
quantization  conditions  to  the  corresponding  wavevector  components  of  the  carriers  To 
model  QCEs  on  the  absorption  edge,  a  photo-excited  electron-hole  two-body  system 
needs  to  be  considered.  In  semiconductors  with  wide  bandwidth,  as  manifested  in  the 
experiments  (see  Fig  12  and  Table  6).,  QCEs  have  a  stronger  j?  blue-shift  effect  from 
the  kinetic  energy  than  the  j-  red-shift  effect  from  the  Coulomb  interaction,  where  L  is 
the  confining  size  of  the  quantum  structure.  As  a  consequence,  QCEs  have  a  stronger 
contribution  from  bandstructure  than  that  of  the  excitonic  effect  despite  the  fact  that 
Coulomb  correlation  modifies  the  optical  spectrum  strongly,  especially  evident  at  low 
temperature.  It  is  essential  to  take  into  consideration  of  the  full  zone  bandstructure  of  the 
semiconductor  when  dealing  with  strong  QCEs  because  approximate  bandstructures,  such 
as  the  Luttinger  model  and  (non-)parabolic  conduction  band  model,  are  far  from  accurate 
even  at  5%  away  from  the  high  symmetry  T  point  (refer  to  Figs.  10  and  1 1).  This  is  the 
very  reason  why  Efros-Efros  theory  [72]  is  deficient  in  treating  strong  QCEs,  as  exhibited 
in  the  experiments  by  Potter  e.t  al.  [58].  Table  4  lists  the  calculated  quantum  dot  sizes 
when  200  meV  discrepancy  between  the  T  transition  points  (refer  to  Eq.  29)  as  predicted 
by  the  TCA  method  occurs  when  using  the  full-zone  bandstructure  and  approximate 
models  for  conduction  and  valence  bands.  Hence,  the  contribution  from  momentum 
quantization  dominates  QCEs  for  ordinary  semiconductors  though  complications  from 
the  excitonic  or  surface  effects  may  be  expected  in  narrow  bandwidth  semiconductors. 
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Table  4:  QD  Sizes  with  200  meV  Discrepancy  between  Calculated 
Absorption  Edges  by  30  x  30  k  ■  p  Model  and  Approximate  Band  Models. 


Material 

Effective  Mass,  m* 

a  (1/eV)  a 

Size,  L  (A) 

GaAs 

0.067 

0.64 

49.6 

CdTe 

0.11 

0.45 

34.8 

ZnS 

0.28 

0.14 

24.4 

ZnSe 

0.14 

0.26 

28.3 

a  is  used  in  the  non-parabolic  band  model  for  the  conduction  band,  as 

h2k2 


E{l  +  aE)  =  - 


2mem* 


The  experimental  results  on  CdTe  quantum  dots  in  a  transparent  glass  matrix  by 
Potter  et  al.  [58]  indicate  that  QCE  on  the  absorption  edge  deviates  from  both  fully 
correlated  and  uncorrelated  excitonic  behaviors  according  to  Efros-Efros  theory  after 
size  distribution  and  finite  barrier  effect  have  been  taken  into  account.  The  blue-shift 
effect  on  the  absorption  edge,  about  1.0  eV  at  the  smallest  particulate  size,  due  to  kinetic 
confinement  is  slightly  offset  by  the  red-shift  effect  due  to  Coulomb  correlation,  which 
can  be  estimated  approximately  using  the  binding  energy  of  a  free  exciton  in  the  bulk 
semiconductor,  i.e., 


&P.T.     


fie 

ot2    o    2 

in  c^Cq 


-  -10.85  meV 


(10) 


with  n  a  mS"mk  and  value  of  parameters  available  in  Table  1.  Actually,  the  binding 
energy  increases  in  lower  dimensions  and  analytical  result  for  2D  yields  a  fourfold 
increase.     Therefore,  the  experimental  data  support  the  simple  argument  of  kinetic 
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confinement  effect  predominating  over  Coulomb  correlation  effect.  As  a  result,  the  TCA 
approach  is  applied  to  CdTe  quantum  dots  system  without  the  Coulomb  effect  to  explain 
QCE  on  the  absorption  edge,  and  a  satisfactory  theoretical  fit  to  the  experiments  is 
obtained  for  T  point  absorption  energies  and  the  fit  for  L  point  energy  E\  is  within  the 
theoretical  error. 

The  EFA  theory  is  applied  to  three  typical  multi-band  Hamiltonian  models  in  order 
to  understand  its  applicability  and  to  compare  it  directly  with  the  TCA  approach.  The 
underlining  issue  is  how  to  interpret  the  results  of  EFA  when  it  is  applicable  to  quantum 
confined  systems.  Concerning  mis,  analytical  solutions  are  sought  for  the  two  generic 
2  x  2  models  under  the  EFA  theory  so  that  ambiguity  can  be  eliminated.  The  effort 
provides  a  satisfactory  explanation  of  the  failure  of  attempt  to  treat  QCEs  by  applying 
the  EFA  theory  with  the  empirical  30  x  30  k  ■  p  model  by  the  author. 

The  present  part  is  organized  as  follows:  (1)  in  Chapter  8,  the  TCA  method  and  the 
EFA  theory  are  described,  and  then  three  models  are  discussed  summarily;  (2)  in  Chapter 
10,  theoretical  results  are  presented  for  both  theories  and  comparisons  are  made  with 
experimental  results;  (3)  in  Chapter  10,  discussions  follow  up  to  make  the  comparison 
of  the  two  theories  more  transparent;  and  (4)  finally,  conclusions  are  given. 
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CHAPTER  8 
THEORIES 

Low  energy  excitations  in  semiconductor  materials  are  usually  satisfactorily  described 
by  Maxwell-Bcltzmann  statistics,  which  results  in  successful  application  of  independent 
electron  quantum  theory  to  quantum  structures.  As  a  note  for  mesoscopic  metallic 
systems,  new  statistical  model  needs  to  be  constructed  because  finite  size  and  coherent 
nature  of  these  systems  invalidate  ensemble  averaging  [62], 

Truncated  Crystal  Approximation 

For  single  electron  in  a  uni-directionally  (z)  confining  semiconductor  structure,  the 
following  Schrodinger  equation  is  solved 

H<f>(r±,z)  =  E<f>{f±tz),  (11) 

where 


with  confinement  potential 


*=&  +  V&  +  7-^-/vV(r-)  x  P 


(12) 


V(f)  '---  {  Vmtl7(f)  ,  0  <  z  <   L  .  (13) 

Vmtn{r)  ,  z  <   L 

Potentials  Vmtli(r)  ,  i  =  1,2,3  are  periodic  in  the  perpendicular  plane  and  assumed  not 
to  undergo  structural  change  in  the  confining  direction  (z).  Under  the  TCA,  the  quantum 
structure  with  potential 

T//^.       ( oo  ,      0  >  z  >  L 
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is  investigated,  which  represents  impenetrable  barriers. 

To  construct  a  complete  TCA  basis  set,  an  observation  is  needed,  which  is,  within 
the  confinement  region,  the  Schrodinger  equation  is  still  satisfied  by  the  solutions  for  the 
bulk  system.  So  the  solutions  for  the  confined  system  are  linear  superposition  of  those 
degenerate  solutions  for  the  bulk  system  and  the  combination  depend  upon  the  broken 
symmetry  from  the  bulk  system  to  the  confined  system.  Therefore,  for  semiconductors  with 
zincblende  structure  (cubic  symmetry),  the  complete  TCA  basis  set  can  be  constructed, 
within  a  normalization  factor  for  uni-directional  confinement  in  the  z  direction,  as 

TCA(r1=f^J_+^(O-^j__^(r-0      ,0<z<L 
*«*JJV  '     \0  ,Q>z>L  {     ' 

where  i>n%{r)  is  the  Bloch  wavefunctions  and  kz  =  kzz  is  given  by 

k*-jj         ,j  =  1,2,  ...,JL  (16) 

with  jL  as  the  maximum  integer  so  that  k  =  k±  +  kx  is  still  inside  the  first  Brillouin 
Zone  (Brillouin  Zone  is  used  strictly  in  place  of  the  first  Brillouin  Zone  throughout  the 
present  thesis). 

The  TCA  method,  in  the  lowest-order  approximation,  assumes  that  solutions  for  the 
wavefunction  <f>(f_i,z)  are  simply  given  by 

<Mf±,*)  =  'P™£(f±,z)  (17) 

so  energy  eigenvalues  are  given  as 

E  =  En(k±  +  jjz^,        j  =  l,2,---,jL.  (18) 

Note  that  at  kj_  =  0,  owing  to  both  quantization  of  the  wavevector  component  kz  and 
nodal  planes  of  the  periodic  part  of  a  Bloch  wavefunction  at  the  confining  barrier[68], 
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the  boundary'  conditions  are  satisfied.  This  basis  set  only  approximately  satisfies  the 
boundary  conditions  away  from  k±  =  0  because  of  the  inequivalency  of  the  two  periodic 
parts  of  the  Bloch  wavefunctions,  «n£x+^  and  Unk±-L-  This  is  one  waY  in  whicn  e,Tors 
are  introduced  into  the  TCA  method.  Using  the  complete  basis  set,  Eq.  (15),  better 
approximations  to  solutions  for  the  Schrbdinger  equation,  Eq.  (11),  with  confinement 
potential,  Eq.  (14),  may  be  constructed.  The  lowest-order  TCA  method  as  outlined 
here  can  be  conveniently  extended  to  three  dimensional  quantum  confinement  cases  for 
semiconductors  with  cubic  symmetry.  Therefore,  the  TCA  method  obtains  quantized 
energies  simply  by  replacing  the  wavevectors  in  the  confinement  directions  into  the  bulk 
bandstructures  with  their  corresponding  quantized  values,  and  then  continuation  needs  to 
be  considered  to  map  out  those  subbands  that  were  previously  not  fully  inside  the  first 
BZ  in  3D  onto  the  new  Brillouin  Zone  in  lower  dimensions  in  order  to  relate  the  total 
DOS  of  ail  subbands  derived  from  each  3D  band  to  the  total  number  of  unit  cells  in 
the  confined  system. 

The  TCA  solutions,  Eq.  (17),  may  be  used,  when  correlation  effects  are  included,  in 
calculations  of  interband  and  inter-subband  optical  transition,  which  are  used  for  detailed 
comparison  with  optical  measurements  of  quantum  confined  systems  where  excitonic 
effects  strongly  modify  optical  characteristics.  Also,  it  is  worth  mentioning  that  a  more 
detailed  analysis  is  needed  to  study  properties  related  to  valence  band  mixing  when 
crossover  between  subbands  occur.  In  the  case  of  multiply  degenerate  satellite  valleys 
in  indirect  bandgap  semiconductors,  the  TCA  method  can  be  constructed  at  local  band 
minima  but  coupling  between  those  accidentally  degenerate  states  within  each  valley  have 
to  be  taken  account  in  the  lowest  order  approximation. 
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Envelop  Function  Approximation 

In  dealing  with  physical  systems  that  are  characterized  by  two  distinctively  different 
length  or  time  scales,  adiabatic  approximations  can  be  used  for  analytical  predictions.  The 
Wenzel-Kramers-Brillouin  (WKB)  theory  in  quantum  mechanics  is  a  good  example  and 
the  effective  mass  approximation  (EMA),  which  is  literally  the  EFA  applied  to  the  lowest 
conduction  band  of  a  semiconductor,  is  another.  However,  adiabatic  approximations  are 
not  supposed  to  work  for  systems  with  comparable  physical  scales,  such  as  for  strong 
confinement  systems  in  which  unit  cell  scale  and  confinement  scale  are  at  the  same 
order  of  magnitude,  i.e.,  ~  10  A.  This  is  precisely  the  reason  why  the  envelop  function 
approximation  fails  in  strong  confinement  regime. 

The  envelop  function  approximation  (EFA),  or  the  approximation  of  the  wave  envelop 
function  [70],  in  the  literature  mainly  treats  conduction  and  valence  bands  independently 
[14,  73]  when  dealing  with  semiconductor  heterostructures  and  superlattices.  However, 
as  suggested  in  White  and  Sham's  work  [74],  modeling  conduction  and  valence  bands  on 
the  same  footing  is  essential  for  better  matchup  of  boundary  conditions,  which,  in  return, 
means  better  results.  Furthermore,  with  the  increase  of  QCE,  dipole  matrix  elements 
between  highest  occupied  states  and  lowest  unfilled  states  can  not  be  predicted  reliably 
by  independent  treatment  of  conduction  and  valence  bands.  Therefore,  when  formulating 
the  EFA  theory  and  selecting  models  in  this  research,  conduction  and  valence  bands  are 
treated  within  the  same  model.  Also,  to  a  certain  limit,  numerical  improvements  are 
expected  when  more  basis  functions  or  bands  are  considered  in  the  model. 

Now  consider  an  electron  moving  in  a  quantum  semiconductor  structure  according 
to  the  Schrodinger  equation,  Eq.  (11),  under  the  potential  given  in  Eq.  (13).  The  first 
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assumption  in  the  EFA  theory  is  that  semiconductors  involved  have  the  same  periodic 
basis  wavefunctions  unjj(r)  {n  =  1,2, •••,00}  at  the  Brillouin  Zone  center1,  and  the 
second  assumption  is  that  they  have  the  same  crystal  structure  in  the  perpendicular  plane. 
The  first  assumption  is  applicable  only  if  no  charge  transfer  occurs  in  the  system,  which 
then  needs  self-consistent  treatment  of  a  correlated  electronic  subsystem,  while  the  second 
one  usually  is  satisfied  although  strain  effects  may  be  introduced. 
The  basis  wavefunctions  unfi{r)  satisfy  the  following  equation 

[P2  +  Y&W  -  «*%!«  -P.        i-  1,2,3  (19) 

for  all  three  materials  in  the  whole  real  space,  and 

WO  -  V&M  +  V3TW,  i-«  1,2,3  (20) 

where  V^j^f)  is  the  antisymmetric  part  of  the  lattice  potential  Vmtn(r)  with  respect  to 
exchange  of  anion  and  cation  so  it  vanishes  for  semiconductors  of  Diamond  structure. 
Using  a  complete  basis  set  of  the  Bloch  wavefunctions 

{"7M         '     *€*•*•.  n-1.*.-.*)},  (21) 

a  single  electron  wavefunction  can  be  expanded  as 

<fi(r±,z)  =  tik?Y,Xn{z)un^f)  ,  (22) 

n 

where  Xn(z)  is  called  the  envelop  function.  Since  the  perpendicular  wavevector  is  still 
a  good  quantum  number,  the  Schrodinger  equation  can  be  reduced  to  its  final  form, 


1.  As  a  new  development  in  the  EFA  theory.  Foreman  [Phys.    Rev.    B54,   1909-21,  1996] 
formulated  the  theory  "with  material-dependent  basis  functions". 
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Eq.  (24),  after  integrating  out  the  fast-aherr.ating  periodic  potential  using  the  following 
approximation 

fdff(f)a(r}^jdff(f}-^      f   drct{f)  ,  (23) 

■unit  cell 
where  function  a(r)  has  the  periodicity  of  the  potential,  and  function  /(f)  is  slowly 

varying  on  the  length  scale  of  function  a(f).  For  semiconductor  heterostructures,  though 

abrupt  changes  in  the  potential  at  the  interfaces  occur,  it  is  often  assumed  that  the  above 

approximation  holds.  The  validity  of  the  EFA  is  called  into  question  with  an  increase 

of  the  degree  of  confinement.    And  the  final  form  of  the  Schrodinger  equation  under 

the  EFA  is 

(p:  +  Pl  +  em-E)Xrn(z)  +  Y{(m\Vasym(r)  +  HS0\n) 

(24) 

+  2{m\P±\n)  ■  PL  +  2{m\Pz\n)  ■  Pz}Xn{z)  =  0  , 

where  (m\0\n}  =       /      dfu*  Jf)Oun^{f)  for  arbitrary  operator,  and  Hso  is  the  Spin- 
unit  cell 

orbital  coupling  interaction.  Note  that  Eq.  (24)  holds  in  different  regions  of  the  system 
with  different  material  parameters  and  is  applicable  to  all  quantum  structures  that  are 
not  under  strong  confinement.  Therefore,  the  envelop  function  approximation  means 
replacing  the  symmetry -broken  quantum  number  with  its  corresponding  operator  in  the 
model  Hamiltonian  of  matrix  form. 

Hamiltonian  Models 

Under  the  envelop  function  approximation,  multiband  models,  which  use  a  limited 
number  of  basis  functions,  may  allow  nontrivial  results,  in  reference  to  the  inverse  square 
behavior  as  predicted  with  a  non-parabolic  model  for  conduction  band  and  Luttinger 
model  for  valence  bands.  Multiband  models  that  treat  conduction  and  valence  bands  on 
the  same  footing  are  discussed  here,  as  mentioned  in  the  previous  section. 


50 


CO^C\JOC\2^'COCOOC\2 

I  I  I  I  *H         *H 

I  I 


(/\9)  A'§J9U3 


■s 

i 

to, 

o 


I 


3 

I 

O 

o 

i 

I 

I 

c 
'n 


3 

4— * 

o 

3 


OQ 

s 


51 


Table  5:  k  ■  p  Results  for  CdTe  Compared  with  Reported  Values  in  Literature. 


Parameter 


Eg  (eV) 


A0  (eV) 


Ai  (eV) 


EvL  (eV) 


Evx  (eV) 


ml 


111  >*mt 


hh 


<  100  >b  ml 


hh 


<  111  >a  ml 


Ih 


<  100  >b  m 


lh 


<lll>«mi 


<  100  >b  m* 


k  ■  p  Results 


i  .603 


0.813 


0.42 


1.189 


2.082 


0.141 


1.258 


0.646 


0.167 


0.190 


0.317 


0.330 


Literature  [29] 


1.59 


0.81 


a    iim  \h2/d2E/dk2\,  defined  for  k  along  T  and  L. 
b  Urn  \h2/d2E/dk2\,  defined  for  k  along  f  and  X. 


fc-o 


0.6 


0. 


3D  30  x  30  k-p  Theory 


Cardona's  15  x  1-5  k  •  p  model  [75-78]  is  extended  to  the  current  30  x  30  one 
by  incorporating  the  spin-orbital  coupling  interaction.  Due  to  a  lack  of  translational 
symmetry,  all  k-p  models  without  modification  are  intrinsically  incapable  of  dealing  with 
QCEs  under  the  EFA  theory.  Therefore,  they  are  restricted  in  use  to  the  TCA  method. 
However,  they  work  satisfactorily  with  the  TCA  method  even  when  strong  quantum 
confinement  is  in  effect,  as  the  results  of  this  research  will  demonstrate  Reconstructed 
k  ■  p  models  that  are  capable  of  dealing  with  quantum  heterostructures  are  available  [79, 
80]. 
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Listed  in  Tables  3  and  5  are  the  complete  set  of  30  x  30  k  ■  p  parameters  and  an 
incomplete  list  of  bandstructure  and  material  parameters  obtained  from  the  computation. 
Figure  13  shows  the  electronic  structure  of  bulk  zincblende  CdTe  as  computed  by  the 
empirical  30  x  30  k  ■  p  method.  It  is  necessary  to  stress  that  all  computation  is  done  at 
zero  temperature,  and  this  could  causes  deviation  from  room  temperature  values  as  much 
as  0.2  eV  as  far  as  energy  levels  are  concerned. 
ID  2  x  2  k  ■  p  Model 

A  2  x  2  k  ■  p  model  is  considered  here.  It  is  analytically  tractable  under  the  EFA 
theory.  This  model  is  used  to  provide  insight  into  the  failure  of  the  30  x  30  k  •  p  model 
under  the  EFA  theory.  It  elucidates  features  that  are  not  transparent  in  the  30  x  30  model, 
such  as  the  origin  of  the  'spurious  states'  and  model-independent  properties.  Furthermore 
these  two  models  allow  direct  comparisons  between  the  EFA  theory  and  TCA  method. 

The  Hamiltonian  for  ID  2x2  fc-  p  model  reads,  in  atomic  units,  as 

ff  =  <i>*    *  +  *»)    •  <25> 

with  Eg  and  P  as  bandgap  and  k  ■  p  coupling  constant,  and  given  as  in  the  3D  30  x  30 
k  ■  p  model.    Eigenvalues  for  Eq.    (25)  are 


*±  =  f  +  *2±j(f)Vw     .  (26) 


ID  White-Sham  Model 


White-Sham  model  was  introduced  for  heterostructures  in  order  to  incorporate  evanes- 
cent states  [74],  or  "wing  states'  to  better  satisfy  interfacial  continuity  conditions.  Derived 
from  the  original  three-band  Kane's  model  [81,  82],  the  current  model  considers  a  two- 
band  effective  Hamiltonian  for  electron  and  heavy  hole,  and  it  is  called  the  'White-Sham 
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model'.  This  model  is  complementary  to  the  previous  model,  but  differs  from  the  previ- 
ous one  in  modeling  the  valence  band  with  a  negative  curvature  (see  the  insets  in  Figs. 
18  and  20  for  details).  The  model  reads,  in  atomic  units,  as 

-k2         Pk     , 


#  = 


[Pk     Ea  +  k2>  W 


with  parameters  Eg  and  P  of  same  meaning  as  in  2  x  2  k  •  p  model.   Eigenvalues  for 
Eq.     (27)  are 

e±  =  f ±  ]j  ( f ) " +  [E* +  p2  +  p)k2  ■  ™ 


CHAPTER  9 
RESULTS 


Oui  theoretical  calculations  include  various  quantum  confinement  structures,  espe- 
cially quantum  dots.  The  results  are  shown  in  Fig.  14-17  with  discussions  followed  in 
next  chapter.  For  comparison,  Table  6  lists  the  relevant  experimental  data  as  appeared 
in  [58]. 


Table  6:  Results  of  CdTe  QDs  and  Calibrated  Absorption  Edges  by  Potter  [58]. 


Average 
Diameter  (A) 

Measured 

F-point 

Absorption 

Edge  Energy 

(eV) 

Calculated 

T-point 

Absorption 

Edge  Energy 

(eV) 

Measured 
L-point  Ei 

Energy  (eV) 

Measured 

L-point  Ei+A] 

Energy  (eV) 

24 

2.05 

2.50 

3.55 

4.14 

30 

203 

2.24 

3.52 

4.16 

43 

1.87 

1.88 

3.43 

4.03 

46 

1.76 

1.79 

3.39 

4  00 

66 

1:68 

1.66 

3.34 

3  92 

79 

1.65 

1.63 

85 

1.63 

1.62 

3.36 

3.97 

86 

1.61 

1.62 

3.32 

3.9 

125 

L56 

1.56 

3.29 

3.88 
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Truncated  Crystal  Approximation 

As  pointed  out  in  the  theory  section,  the  TCA  method  for  regular  quantum  structures 
means  mapping  3D  electronic  structural  cross  sections  into  a  lower  dimensional  BZ  along 
the  quantization  direction(s).  Results  of  ID  quantum  confinement  effects  on  the  electronic 
energy  levels  for  CdTe  quantum  film  are  shown  in  Fig.  14.  As  clearly  seen,  the  warping 
nature  of  the  valence  bandstructure  results  in  negative  hole  effective  masses  for  certain 
subbands  when  the  confinement  effect  is  strong. 

In  order  to  model  realistic  (spherical)  quantum  dots  which  are  confined  in  3D,  an 
equal  volume  scheme  is  adopted.  Calculation  of  bandstructure  quantization  of  a  quantum 
dot  is  done  for  a  quantum  cube  of  the  same  volume  of  the  spherical  dot  as  measured 
in  the  experiments  [58].  That  is,  the  confinement  width  Lcube  is  scaled  according  to 
Lcvbe  =  i/^D^t  with  Ddot  being  the  diameter  of  the  quantum  dot.  This  scheme  works 
in  the  loose  sense  of  interpolation.  Then,  for  the  sake  of  continuation  because  there  is 
no  Brillouin  zone  anymore,  association  of  quantized  wavevectors  with  high  symmetrical 
points  in  the  Brillouin  Zone  is  made  as  follows, 

,7T      7T      7T 

<ZTT>         -  r 
iVVn*T)       -X 

,*    .       *  x  (29) 

{T3l*Tjl*l]  "*  K 

UL  *  j,]L  *  j,JL  *T)~*    L  ' 

where  ji  is  the  maximum  integer  allowed  so  that  the  corresponding  wavevector  is  within 
the  Brillouin  Zone.  For  different  symmetric  points,  ji  differs. 

Comparison  of  the  theoretical  results  with  the  experimental  optical  transition  energies 
of  CdTe  QDs  in  a  transparent  glass  matrix  [58]  is  shown  in  Fig.    17.   Over  the  whole 
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Figure  14:  Truncated  Crystal  Approximation  (TCA)  result  of  30  A  CdTe  quantum  film. 

The  TCA  method  maps  the  quantized  cross  sections  of  the  electronic  structure 

onto  the  unaffected  plane.  Shown  are  the  first  three  subbands  from  each 

band  in  the  bulk.  Reflected  clearly  are  the  warping  nature  of  the  valence 

bands  and  the  lift  of  the  double  degeneracy  at  the  top  of  the  valence  bands. 
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QUANTUM   CONFINEMENT 
TCA  RESULTS   for   CdTe 
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Figure  15:  Predicted  bandedges  at  symmetrical  points  of  CdTe  quantum  cube  by 
the  TCA  method.  No  inverse  square  behaviors  are  found  and  asymptotic 
saturation  of  quantum  confinement  effect  on  the  bandedges  is  predicted. 
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Figure  16:  Predicted  transition  energies  of  CdTe  quantum  cube  by 
the  TCA  method.  V  point  is  the  absorption  edge  above  16  A  and 
then  replaced  by  an  unidentified  point  in  some  general  direction. 
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range  of  measurement,  the  resuh  ,  ,„  e*ce,ie„,  „„,  wWl  lhe  cxpenmenI    ^ 
experiment.,  data  has  been  corrected  fof  ^  o|  _  ^^  ^  ^  ^^ 

into  glass  barriers. 

Despite  the  face  that  the  TCA  method  is  sat.sfaconiy  applied  t„  [he  CdTe  ouantum 
confines  problem,  the  appHcat,on  of  (1k  EpA  me(hod  wth  m  ^  m  .  _  ^  ^  ^ 

probiern  turns  out  ,„  yield  „„.physlcal  resu„,   To  ^^  ^  ^^  ^  ^  ^  ? 
analydca.  Ham,„on,an  models  are  used  Wlth  bM|,  methods  for  ^^   ^  ^ 
are  show,,  in  the  next  section  together  with  those  of  the  EFA  method. 
The  TCA  solutions  for  ID  2  x  2  i  ■  p  mode|  are 

As  shown  ,n  Figs.   ,8  and  ,9,  no  „„,  of  any  klnd  „  produced  ^  fc  ^ 

method  a„d  TCA  approach.  However,  for  the  White-Shan,  mode,,  whose  TCA  so.utions 
are 


e 


l-h0) +*.+»+ 


Excellent  agreement,  as  shown  in  Fig    20   is  fn.mH  » 

rig.    ^o,  1S  found  among  analytical,  numerical,  and 


TCA  results. 


Envelop  Function  Approximation 

to  eincdate  the  failure  of  the  appi.cation  of  ,he  EFA  method  with  30  x  30  f  •  p  mode, 
The  EMA  theory,  the  simpiest  form  of  the  EFA  theory,  predicts  the  we^nown  inverse 
soaare  dependence  of  quantization  of  energy  ,evels  npon  confinement  width.  The  same 
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Figure  17:  Experimental  data  compared  with  TCA  predictions  for 

CdTe  quantum  dots.  The  experimental  data  are  corrected  by  considering 

size  distribution  and  tunneling  into  finite  glass  barriers.  Discrepancy  at  large 

L  limit  is  because  the  TCA  predictions  are  for  QDs  at  zero  temperature. 
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inverse  square  behavior  is  found  for  the  White-Sham  model,  as  shown  below.  Though 
improved  in  agreement  through  more  accurate  bandstructures,  the  2  X  2  k  ■  p  model  does 
not  yield  converging  results. 

Parameters  for  the  2x2  models  are  fit  to  zincblende  ZnS  instead  of  CdTe  but  con- 
clusions are  not  affected.  Finite  difference  scheme  is  adopted  to  solve  the  corresponding 
differential  equations  for  the  2  x  2  models.  Numerical  and  analytical  results  are  obtained 
simultaneously  for  assurance. 

The  analytical  results  for  2x2  k  •  p  model  are  obtained  by  solving  Eq.  (32) 
numerically  for  the  quantized  wavevectors  ku  k2,  and  the  quantized  energies  t.  The 
algebraic  equations 


\smkiLs'mk2L  -  (1  -  cos  kiL  cos  koL)*  (32) 

2Al(k1)Al(k2)A2  (k1)A2(k2 ) 

*I(*i)A3(te)+A]!(*a)&|(]k, ) 
where  Ax(fc)  =  Eg  +  k2  -  e,  A2{k)  =  -Pk,  for  quantized  energies  t  at  k}  and  k2  are 
derived  when  solving  an  eigenvalue  problem  under  zero  boundary  conditions  [63].  The 
analytical  results  are  compared  with  TCA  results  in  Figs.  18  and  19,  and  no  agreement 
between  the  two  methods  is  obtained.  Also  diverging  behavior  is  found  for  this  model 
when  numerical  solutions  are  sought.  This  behavior  will  be  addressed  in  the  discussion, 
though  it  is  worth  pointing  out  that  the  valence  band,  in  the  whole  real  plane,  introduces 
'spurious  states',  i.e.,  those  larger  k  solutions  in  Eq.  (32)  when  t  >  Eg,  which  are 
unphysical  and  due  to  the  faulty  model.  Because  of  the  lack  of  numerical  convergence 
in  2  x  2  k  •  p  model,  the  effort  to  extend  the  Envelop  Function  Approximation  to  30  x  30 
k  ■  p  model  was  abandoned.  In  retrospective,  the  study  of  the  2  x  2  %  •  p  model  was  to 
shed  light  on  the  cause  for  the  numerically  diverging  behavior. 


62 


ID    2x2   k.p   model 


6.0 


5.0     - 


v 


4.0      - 


mesh   points^    4 
P=1.2,    Eg  =  .265 


3.0 


b  2° 

tic 

03 

G      1.0 


0.0 


-1.0 


.  i-j-ft  r~f-*~~'  *'TTtTtTtT i..i- 

45.0 
^35.0 
%  2.5.0 
""15.0 

w    5.0 
-5.0 


0.00         0.10 
k   (    27r/ao 


0.0     0.3     0.6     0.9     1.2     1.6     1.9     2.2     2.5 
WIDTH"2       (1/A"2)    (x    10"3) 


Figure  18:  Comparison  between  numerical  EFA  and  TCA  results.  Solid  lines 

are  EFA  calculations  with  4  mesh  points  in  finite  difference  scheme.  Inset 
shows  ID  2  x  2  k   p  model,  which  is  asymptotically  parabolic  for  both  bands. 
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Figure  19:  Comparison  between  numerical  EFA  and  TCA  results.  Solid  lines  are  for 

EFA  calculations  with  20  mesh  points.  Compared  to  EFA  results  with  4  mesh 

points  (Fig.   18)  and  TCA  results,  no  agreements  between  EFA  results  with 

different  mesh  points  and  between  EFA  and  TCA  results  have  been  achieved. 
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The  White-Sham  mode!  is  studied  in  comparison  with  the  same  parameters  for  2  x  2 
k  ■  p  model  in  order  to  understand  the  differences  among  various  models.  The  quantized 
energy  and  the  wavevectors  k\,  k>  are  obtained  by  solving  the  following  algebraic 
equations 

sin  kiLs'm  koL  =  (1  —  cos  k\L  cos  k->L)*  @3) 

2A1(fc1)A1(fc2)A2(fe1)A2(fc2) 
AJ(Jki  )^(ij)+ AJ(*»)  A*  (*, ) 

where  Ai(fc)  =  £9  +  fc2  -  e,  A2(fc)  =  -PA;.  The  difference  between  2  x  2  k  •  p  model 
and  the  White-Sham  model  is  the  way  in  which  the  valence  band  is  modeled:  in  the 
former  model,  the  valence  band  is  an  asymptotical  positive  parabolic  which  allows  real 
solutions  to  the  other  pair  of  fc's  in  Eq.  (32)  when  e  >  Eg;  while  in  the  White-Sham 
model,  the  the  valence  band  is  an  asymptotical  negative  parabolic  which  allows  only 
imaginary  solutions  to  the  other  pair  of  fc's  in  Eq.  (33)  when  e  >  Eg.  This  distinctive 
difference  is  obvious  in  the  two  insets  of  Fig.  18  and  Fig.  20.  As  shown  in  Fig.  20, 
excellent  agreement  is  achieved  between  the  EFA  and  TCA  method. 
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Figure  20:  Analytical,  numerical  and  TCA  results  of  the  White-Sham  model.  Straight 
solid  lines  are  for  viewing  guide.  The  inset  shows  the  bands  as  predicted  by  the 
White-Sham  model.  Excellent  agreements  have  been  achieved  among  all  results. 


CHAPTER  10 
DISCUSSION 

Multiband  models  using  the  EFA  theory  simply  mix  degenerate  states  from  the  bulk 
material  when  dealing  with  QCEs,  and  a  mixture  state  in  quantum  structure  is  far  from 
similar  to  an  envelop  function  on  top  of  the  original  Bloch  wavefunction.  Hence  it  is 
reasonable  to  believe  that  the  EFA  theory  is  at  fault  in  dealing  with  QCE  because  of 
this  lack  of  self-consistency.  But  every  evidence  shows  that  this  belief  is  far  from 
the  truth.  This  is  an  example  of  a  non-self-consistent  theory  that  works.  This,  in 
another  way,  suggests  why  the  TCA  method  succeeds  when  applied  to  QCE  problems 
because  those  basis  functions,  Eq.  (15),  constructed  based  upon  the  notion  of  destructive 
interference,  approximately  lake  into  consideration  both  degenerate  states  mixing  and 
boundary  conditions. 

Truncated  Crystal  Approximation 

When  dealing  with  quantum  confinement  systems,  because  of  the  abrupt  changes 
in  the  potential  profiles,  the  EFA  theory  will  not  only  couple  degenerate  states  but 
also  discretize  the  electronic  energy  levels  in  order  for  the  wavefunctions  to  satisfy 
the  boundary  conditions.  In  order  to  handle  the  ever-decreasing-in-size  nanosystems 
empirically,  more  and  more  basis  functions  are  needed  to  construct  an  accurate  multiband 
Hamiltonian.  It  means  more  difficulty,  if  not  impossible,  to  solve  the  coupled  differential 
equations  either  analytically  or  numerically  under  the  EFA  method. 

On  the  other  hand,  though  unfit  for  treating  quantum  confinement  problems  under 
the  EFA  theory  because  of  'spurious  states',  30  x  30  k-p  model  under  the  TCA  approach 
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has  been  successfully  applied  to  treat  quantum  confinement  problems.  Surprisingly  as  it 
seems,  the  reason  for  this  paradox  is  as  follows.  All  models  that  have  been  discussed 
in  this  research  are  approximate  ones  since  all  of  them  fail  to  include  the  fundamental 
symmetry  of  spatial  periodicity  in  solids.  The  Hamiltonians  would  have  been  periodic  in 
k  space  if  translational  symmetry  were  considered.  This  very  defect  in  the  k  ■  p  models 
introduces  'spurious  states',  which  are  unphysical  solutions  and  beyond  the  Brillouin 
Zone,  when  QCEs  are  treated  under  the  EFA  theory.  It  is  this  artificial  deficiency 
in  the  k  ■  p  model  that  spells  failure  for  application  of  the  EFA  theory  to  QCE.  For 
the  TCA  approach,  it  is  the  very  approximation  that  was  made  when  constructing  the 
basis  functions  to  decouple  states  in  the  Brillouin  Zone  from  those  degenerate  states 
beyond  the  Brillouin  Zone  fundamentally  guarantees  the  success  of  the  approach.  The 
physical  reason  for  the  success  is  because,  as  pointed  out  in  discussion  above,  the  TCA 
basis  functions,  Eq.  (15),  constructed  upon  notion  of  destructive  interference,  take  both 
degenerate  state  mixing  and  boundary  conditions  into  consideration  [68].  However,  when 
sizes  of  quantum  structures  are  reduced  further,  their  physical  properties  deviates  more 
and  more  from  those  of  bulk  materials  and  approach  those  of  molecules.  This  extreme 
limit  invalidates  the  TCA  method,  which  predicts  the  properties  of  a  cluster  from  those 
of  bulk  material. 

Improvement  of  the  30  x  30  k-p  model  under  the  TCA  method  calculation  is  possible 
if  more  TCA  basis  functions  are  used  in  a  higher  order  approximation.  Furthermore, 
as  suggested  in  Fig.  17,  inclusion  of  excitonic  effect,  which  should  reduce  transition 
energies,  should  improve  the  agreement  between  the  theory  and  the  experiments  as  well. 
It  is  beyond  this  crude  approximation  of  the  TCA  method  to  study  quantum  confined 
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systems  which  change  drastically  and  qualitatively  with  the  size  and  dimension  of  the 
systems.    Then,  higher  order  approximations  are  required,  such  as  in  the  case  of  Si 
quantum  wires  and  thin  films,  which  involves  multiple  degeneracy. 

Envelop  Function  Approximation 

The  experiments  [58]  show  deviation  from  the  j?  inverse  square  behavior,  i.e., 
energy  levels  quantize  in  proportion  to  the  inverse  square  of  the  confinement  size  of 
the  structure.  As  calculation  indicates,  the  rather  general  models  that  are  discussed 
under  the  EFA  theory  either  reproduce  trivial  inverse  square  behavior  or  fail  to  yield 
physically  meaningful  results  over  the  whole  confinement  range.  Beyond  the  weak 
quantum  confinement  range,  first  the  theory  becomes  incapable,  then  the  theory  turns 
invalid  when  the  confinement  size  is  comparable  to  the  unit  cell.  Therefore,  the  envelop 
function  approximation  may  be  applicable  to  weak  quantum  confinement  problems,  but 
care  in  selecting  Hamiltoman  models  should  be  exercised.  Also,  the  Efros-Efros  theory 
predicts  the  same  trivial  inverse  square  behavior  even  though  it  considers  Coulomb 
correlation,  which  indicates  the  inadequacy  of  this  theory  in  dealing  with  QCE  and  the 
less  importance  of  the  excitonic  effect. 

Defect  in  the  models,  which  means  lack  of  consideration  of  the  spatial  periodicity, 
and  the  approximate  nature  of  the  EFA  theory  are  the  reasons  for  the  inadequacy  of  the 
theory  to  treat  quantum  confinement  effect.  For  the  ID  2  x  2  k  ■  p  model,  analytically 
there  are  always  two  positive  solutions  for  k2  according  to  Eq.  (32)  when  e  >  Eg,  so  a 
pair  of  real  solutions  for  k  exist  from  the  larger  solution  of  k2.  This  pair  of  real  solutions 
for  k  are  called  'spurious  states'  [70]  because  they  are  outside  the  Brillouin  Zone  and 
unphysical  since  they  would  be  nonexistent  if  spatial  periodicity  of  the  system  were 


69 
considered.  Numerically,  these  "spurious  states'  are  mapped  as  energy  eigenstates  from 
outside  the  Brillouin  Zone,  thus  cause  numerical  instability  and  convergence  problem. 

Different  from  ID  2  x  2  k  ■  p  model,  for  ID  White-Sham  model,  there  are  always 
one  positive  and  one  negative  solution  for  k2  according  to  Eq.  (33)  when  t  >  Eg.  Then 
the  negative  solution  for  k2  yields  a  pair  of  imaginary  solutions  for  k  which  are  called 
'wing  states'  [74]  because  they  have  larger  amplitudes  than  the  real  solutions  of  k.  These 
'wing  states'  are  evanescent  states,  so  they  cause  no  numerical  problems.  The  two  bands 
in  the  White-Sham  model  asymptotically  decouple  into  two  independent  bands  except 
for  introducing  the  'wing  states',  which  is  equivalent  to  two  independent  bands  under 
the  effective  mass  approximation.  It  is  no  wonder  that  this  model  reproduces  the  trivial 
inverse  square  behavior.  Therefore,  the  White-Sham  model  under  the  EFA  theory  is 
applicable  to  QCEs  but  reproduces  trivial  results  in  strong  confinement  regime,  and  2  x  2 
k  ■  p  model  can  not  be  used  under  the  EFA  theory  due  to  the  defect  in  the  model. 

As  a  note,  the  Luttinger  model  for  valence  bands  [83]  is  also  investigated  and 
the  same  trivial  inverse  square  behavior  for  small  quantum  systems  is  found  though 
better  quantitative  agreement  is  achieved  beyond  intermediate  quantum  confinement 
range  between  the  EFA  theory  and  the  TCA  method,  though  no  results  are  given 
here.  Conclusion  can  be  drawn  that  multiband  Hamiltonian  model  that  does  not  produce 
'spurious  states'  is  applicable  for  study  of  QCEs  under  the  EFA  theory  and  quantitatively 
improves  the  agreement  as  the  number  of  basis  functions  increases,  though  trivial  results 
are  produced  in  strong  confinement  regime. 


CHAPTER   11 
CONCLUSIONS 

Quantum  confinement  effects  on  the  absorption  edge  in  CdTe  quantum  dots  have 
been  satisfactorily  modeled  by  using  an  empirical  full  zone  30  x  30  k  ■  p  model  with 
the  TCA  method  The  theoretical  results  are  in  agreement  with  experimental  data  for 
zincblende  CdTe  clusters  in  a  transparent  glass  matrix.  However,  the  application  of  this 
model  within  the  EFA  theory  framework  fails  to  produce  converging  numerical  results. 

Study  on  k-p  models  and  White-Sham  model  with  both  the  EFA  theory  and  the  TCA 
approach  allows  exploration  of  detailed  aspects  of  these  models  and  comparison  between 
the  two  theories.  The  k  ■  p  models  with  the  TCA  method  are  applicable  to  make  reliable 
predictions  for  quantum  confinement  effects  until  the  molecular  aspects  of  the  quantum 
systems  dominate.  Furthermore,  increasing  the  number  of  basis  functions  improves  the 
accuracy  of  the  modeling. 

On  the  contrary,  under  the  EFA  theory,  owing  to  the  lack  of  consideration  of  the 
translational  symmetry,  k  ■  p  models  are  inapplicable  for  modeling  QCE.  On  the  other 
hand,  the  White-Sham  model,  though  suitable  for  describing  QCE,  predicts  the  trivial 
inverse  square  behavior  in  strong  quantum  confinement  regime,  which  is  a  manifestation 
of  the  asymptotically  parabolic  dispersions  for  large  wavevectors  in  White-Sham  model 
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PART  IV 
STATIC  SCREENING  AND  MOBILITY  IN  n-ZnSe 


CHAPTER  12 
INTRODUCTION 


Zinc  selenide  is  a  commercial  blue/green  light-emitting  diode  (LED)  material,  and  has 
been  intensively  studied  in  recent  years.  In  quantum-well  (QW)  LED  and  laser  devices, 
ohmic  contacts  need  to  be  made,  as  electrodes,  to  the  active  region  where  carriers  can 
be  injected,  recombine,  and  thus  emit  light.  Heavily  doped  ZnSe  is  commonly  used 
in  the  bandgap  grading  scheme  [12]  to  make  ohmic  contact,  and  the  carrier  mobility 
of  intermediately  to  heavily  doped  semiconductor  is  an  important  reliability  parameter. 
Among  all  major  mechanisms  that  determine  the  mobility,  ionized  impurity  scattering 
becomes  predominant  at  low  temperature  (usually  <  25  K),  and  gets  enhanced  under 
heavy  doping.  Furthermore,  heavy  doping  can  induce  metal-insulator  transitions  in 
semiconductors. 

Heavily  doped  semiconductors  have  been  extensively  studied  for  possible  applications 
in  solar  cells,  tunneling  devices  and  metallization.  In  a  heavily  doped  semiconductor, 
impurity  band  may  be  formed  due  to  spatial  overlap  between  the  electronic  wavefunctions 
of  neighboring  impurities.  Disorder  in  the  system  leads  to  the  existence  of  a  mobility 
edge  in  the  impurity  band.  Above  the  edge,  conduction  is  free-electron-like,  and  below 
the  edge,  conduction  is  possible  by  variable-range  hoppings.  A  semiconductor  becomes 
metallic  when  the  chemical  potential  of  the  system  passes  through  the  mobility  edge  with 
increasing  doping  level.  This  transition  is  well-known  as  the  Anderson  transition  [84-86]. 
In  comparison,  intermediately  doped  semiconductors  are  not  well  understood,  especially 
near  the  transition.  This  research  is  motivated  by  this  fact  and  attempts,  qualitatively,  to 
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explain  the  experiments  by  Park  and  Simmons  [87]. 

In  order  to  understand  a  system  with  wide  range  of  doping,  a  self-consistent  theory 
was  proposed  for  determination  of  the  chemical  potential,  static  screening  effects  and 
ionization  of  impurities.  For  simplicity,  only  single-level  hydrogen-like  donors  and 
acceptors  are  considered  for  impurities.  Static  screening  is  treated  by  a  generalized 
Debye  theory1.  Mobility  due  to  ionized  impurity  scattering  is  obtained  by  application  of 
the  relaxation  time  approximation  to  the  Boltzmann  kinetic  equation  (BKE)2.  The  theory 
is  applicable  to  other  widegap  semiconductors,  though  formation  of  impurity  bands  at 
high  doping  concentration  has  not  been  taken  into  account  because  (1)  the  intermediate 
doping  regime  is  of  major  interest  for  this  research;  (2)  the  impurity  bands  may  disappear 
as  a  result  of  screening;  (3)  the  mobility  is  mainly  determined  by  free  carriers  in  the 
bulk  conduction  band.  Also,  we  neglect  the  bandgap  narrowing  due  to  doping,  which  is 
assumed  to  be  far  smaller  than  the  wide  bandgap  of  the  system  [7].  No  consideration 
of  fluctuations  in  impurity  potentials  and  impurity  conduction  has  been  taken.  Only  a 
uniformly  doped  semiconductor  is  under  consideration. 

It  is  found  that  [1}  static  screening,  as  predicted  by  the  generalized  Debye  theory, 
can  induce  an  abrupt  metal-insulator  (MI)  transition,  with  a  critical  temperature  found  at 
26.5  K,  for  uncompensated  systems  with  donor  concentration  above  101 7cm-3  inclusive 
and  no  such  transition  for  a  highly  compensated  system  (see  Figs.  26  and  33);  (2) 
static  screening,  as  predicted  by  the  generalized.  Debye  theory,  induce  a  smooth  MI 
transition  for  highly  compensated  systems  and  the  critical  donor  concentration  is  found 
at  1017 cm-3;  (3]  static  screening  drastically  affects  the  system  when  donor  concentration 

1.  Sometimes  called  rhe  Debye  Huckel  theory 

2.  Interchangeably  used  for  the  Boltzmann  transport  equation  in  the  present  thesis 
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is  well  below  the  critical  concentration  of  metal-insulator  transition,  given  by  Eq.  (34); 
(4)  donor  level3  is  reduced  due  to  static  screening;  (5)  free  earners  can  be  reduced  by 
an  order  of  magnitude  for  compensated  systems  at  high  donor  concentration  (see  Fig 
34);  (6)  the  generalized  Debye  theory  over-estimates  the  static  screening  effect.  The 
present  theory  was  compared  with  the  experiments  of  Park  and  Simmons,  and  a  partial 
and  qualitative  agreement  is  achieved. 


3.  or  donor  binding  energy,  as  used  throughout  the  present  part. 


CHAPTER   13 
BACKGROUND 


Carrier  mobility  in  a  semiconductor  device  is  an  important  parameter  which  essen- 
tially determines  how  fast  the  device  can  operate.  One  of  the  technological  aspirations 
of  semiconductor  community  is  to  increase  the  carrier  mobility.  Higher  mobility  means 
speedier  devices  and  improvements  in  performance.  Mobility  of  free  earners  is  deter- 
mined by  how  often  and  how  intense  they  scatter  off  irregularities  in  a  semiconductor, 
such  as  defects,  impurities,  interfacial  discontinuities  and  dynamical  excitations  Among 
those  important  scattering  mechanisms  that  limit  the  carrier  mobility  in  widegap  semicon- 
ductors near  equilibrium  are  polar  LO-phonon  and  ionized  impurity  scatterings.  Carrier- 
carrier  and  plasmon  scatterings  become  prominent  for  poorly  screened  systems  with  high 
carrier  concentration.  Doping  not  only  introduces  impurities  but  also  increases  carrier 
concentration,  therefore  it  affects  carrier  dynamics.  Overdoped  semiconductors  may  un- 
dergo Anderson  transition  and  become  metallic,  but  with  lower  carrier  mobility.  The 
critical  donor  concentration  Nf  for  uncompensated  n-type  semiconductors  is  predicted 
by 

A7' 4^0.27  (34) 

with  adH  being  the  hydrogen  radius  h  t^KJrn^m^''  for  electrons,  n  being  the  relative 
dielectric  constant1,  eo  being  the  dielectric  constant  in  vacuum,  m*e  being  the  relative 
electron  effective  mass,  and  me  being  free  electron  mass.  Using  parameters  in  Table  7, 


1.  In  Table  1,  it  is  denoted  as  cq. 
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it  is  found,  for  uncompensated  n-ZnSe,  that 


i!7 -3 


Nf  =  8.389  x  101' cm"4.  (34') 

The  overall  performance  of  semiconductor  devices  is  enhanced  if  operated  at  low 
temperature  because  of  reduction  of  the  number  of  scattering  events  due  to  thermal 
excitations.  The  influence  of  different  kinds  of  scatterings  on  carrier  mobility  depends 
upon  temperature  and  carrier  concentration  differently.  Conventionally  speaking,  ionized 
impurity  scatterings  have  a  weaker  temperature  dependence  than  other  mechanisms  but  a 
stronger  carrier  concentration  dependence.  Hence  at  low  temperature  (for  example,  4.2  K) 
or  high  impurity  concentration  (over  I018cm~3),  ionized  impurity  scattering  determines 
the  carrier  mobility  of  a  semiconductor,  while  polar  LO-phonon  scatterings  compete 
with  ionized  impurity  scatterings  as  the  temperature  rises  or  impurity  concentration 
decreases.  For  elastic  scatterings,  which  ionized  impurity  scatterings  are,  relaxation 
time  approximation  can  be  applied  to  predict  the  earner  mobility  which  is  directly 
proportional  to  the  relaxation  time.  The  relaxation  time  is  defined  with  respect  to  energy, 
which  measures  the  intensity  of  scatterings  at  a  given  energy.  However,  optical  phonon 
scattering  at  medium  temperature  can  not  be  treated  as  an  elastic  process,  and  therefore 
a  relaxation  time  approximation  is  invalid. 

Owing  to  its  strong  dependence  upon  carrier  concentration,  the  conventional  weak 
temperature  dependence  of  ionized  impurity  scatterings  can  be  drastically  changed  in  and 
above  intermediate  temperature  range,  given  the  approximately  exponential  relationship 
between  carrier  concentration  and  temperature.  First  of  all,  temperature  and  impurity 
concentrations  are  independent  variables,  while  carrier  concentrations  depend  upon  tem- 
perature and  impurity  concentrations.  Then,  the  chemical  potential  varies  according  to 
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charge  neutrality  at  different  temperatures  Thirdly,  the  carrier  statistics  may  change  from 
Maxwell-Boltzmann  (MB)  to  Fermi-Dirac  (FD)  for  degenerate  carriers,  and  this  transition 
completely  changes  the  screening  characteristics.  The  critical  electron  concentration  nCT 
depends  upon  temperature  and  is  given  by  nCT  =  NcFi_(0)  ,  where  Arc  and  J^iirj)  are 
defined  in  Eq.  (42).  Finally  and  most  importantly,  impurity  states  are  modified  as  carrier 
concentrations  increases,  mainly  due  to  the  static  screening  effect.  This  complexity  alters 
the  conventional  characterization  oi  ionized  impurity  scatterings  and  affects  the  modeling 
of  mobility  in  doped  semiconductors.  In  this  research,  only  the  contribution  from  ionized 
impurity  scatterings  is  considered. 

This  research  treats  carrier  scatterings  due  to  ionized  impurities  by  taking  into  full 
account  both  the  FD  statistics  and  self-consistent  determination  of  the  chemical  potential 
which,  in  return,  determines  carrier  concentrations  and  their  derivatives  with  respect  to 
the  chemical  potential.  The  ionized  impurity  scatterings  are  treated  with  the  relaxation 
time  approximation,  but  without  using  either  the  Brooks-Herring  [88]  or  the  Conwell 
Weisskopf  [89]  approximation.  The  present  theory  ignores  the  formation  of  impurity 
bands,  as  stated  in  the  previous  chapter,  and  its  influence  on  the  static  screening  may 
be  estimated  as  follows.  First  of  all,  the  density  of  states  (DOS)  due  to  impurities  does 
not.  change  after  metal- insulator  transition.  Instead,  there  is  a  gradual  spread  of  DOS 
with  energy  centered  around  the  single  impurity  value,  which  is  equivalent  to  a  system 
of  distributed  impurity  levels.  At  low  temperature,  the  screening  increases  appreciably 
and  modifies  the  present  results.  Secondly,  conduction  by  impurity  states  is  negligible 
at  high  enough  temperature  compared  to  conduction  by  the  bulk  states  because  of  their 
hopping  nature    Finally,  a  large  fraction  of  the  impurity  states  are  excited  into  the  bulk 
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states  resulting  from  the  reduction  of  donor  levels  due  to  increased  static  screening  at 

high  temperatures. 

The  present  part  of  the  research  is  organized  as  follows:  First  of  all,  the  set  of 
equations  for  determination  of  chemical  potential  in  a  homogeneous  semiconductor  are 
summarized.  Then,  the  expression  for  mobility  due  to  ionized  impurity  scatterings  is 
deduced.  Furthermore,  the  results  from  applying  the  theory  to  n-ZnSe  are  presented  and 
compared  to  the  experimental  data  [87].  Discussion  follows  up  and  finally  conclusions  are 
drawn.  (In  Appendix  C,  exposition  of  the  generalized  Debye  theory  and  its  applications 
are  given  ) 


CHAPTER  14 
THEORIES 

Chemical  Potential 


Chemical  potential  is  determined  by  thermodynamical  equilibrium  of  a  system.  For 
homogeneous  semiconductors,  the  condition  for  the  equilibrium  is  reduced  to  charge- 
neutrality  among  ionized  impurities  and  charged  earners.  A  change  in  the  chemical 
potential,  in  turn,  affects  carrier  concentrations,  dynamical  parameters,  such  as  scattering 
rates,  and  static  parameters,  such  as  dielectric  constant.  Table  7  lists  the  material 
parameters  used  in  the  present  calculation  for  n-ZnSe.  The  bottom  of  the  conduction 
band  is  used  as  the  reference  for  the  chemical  potential,  or  EC{T)  =  0.  Also,  all  other 
material  parameters  are  positive  and  in  their  absolute  values  so  that  the  bandgap  Eg 
directly  appears  in  the  following  formulas. 

First  of  all,  listed  below  are  the  three  equations  governing  charge  equilibriums,  which 
hold  for  acceptors,  donors,  and  total  charge  neutrality. 

Na  =  N~  +  Nl ,  (35) 

Nd  =  Nl  +  Nl ,  (36) 

N-  +  n  =  iV+  +  p ,  (37) 

where  Na,  N~ ,  N£,  Nd,  Nf ,  Nf,  n,  p  are  concentrations  of  acceptors,  ionized 
acceptors,  neutral  acceptors,  donors,  ionized  donors,  neutral  donors,  electrons,  and  holes. 
Equilibrium  equations  (35-37)  are  self-evident.  In  this  research,  compensation  ratio  9 
is  defined  as  9  =  N~/N£. 
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Table  7:  ZnSe  Material  Parameters  Used  in  the 

> 

Calculations 

of  Part  IV. 

Parameter 

Notation 

Value 

Reference 

electron  effective  mass 

m* 

0.17 

[90] 

hole  effective  mass 

"»: 

0.60 

[91] 

energy  gap 

Eg  (eV) 

2.70 

[92] 

relative  dielectric  constant 

K 

9.20 

[93] 

donor  level 

E\  (meV) 

27.3a 

acceptor  level 

El  (meV) 

96.4b 

acceptor  hydrogenic  radius 

a%  (A) 

8.1  r 

donor  hydrogenic  radius 

adH  (A) 

28.6d 

■ 

F°  -- 

=  13.6m*  k"2. 

b 

pO   - 

=  13.6m*, /c-2 

c 

«ff    = 

=  .529/c/mJ. 

d 

4  = 

=  .529/c/m*. 

Then,  statistics  for  all  the  thermodynamicaily  active  parties  in  the  system  are  con- 
sidered. FD  statistics  for  carriers  are  used  in  all  parametric  ranges.  The  expressions  are 
listed  below  and  are  available  in  most  solid  state  physics  textbooks,  such  as  Ashcroft 
and  Mermin  [94], 


-V, 


n  =  NcTi(tj)  , 

2 
2 


.  \e-0{Ea- E3-Ef)  ' 

Nd 


N5  =  - 


l  +  \e-f>(Ed+Ef)> 


(38) 
(39) 
(40) 

(41) 
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where 


Nc  =  J™>MT\*  =  4g29  x  1015(Tm:)|cm-3  j 

3 


(42) 


X 


1 


2      /" 


V??  /       1  +  exp  (e  -  77)  ' 

o 

T  being  the  temperature  and  Ef  being  the  Fermi  energy. 

In  combination  with  Eqs.  (35-41),  the  binding  energy  of  a  Hydrogen-like  impurity 
is  treated  by  including  static  screening  using  the  generalized  Debye  theory,  which  is 
derived  in  Appendix  C.  To  self-consistently  solve  the  equations,  the  following  procedure 
is  performed.  First,  of  all,  the  Fermi  energy  Ef  is  determined  by  solving  Eqs.  (35-41) 
with  the  material  parameters  given  in  Table  7  and  variables,  i.e.,  temperature  T,  donor 
concentration  N4,  and  acceptor  concentration  Na,  given.  Then,  the  Debye  screening 
length  \q  is  determined  from 

The  new  binding  energies,  Ed  for  donors  and  Ea  for  acceptors,  are  updated  with  the 
following  relations  [95] 

EA  A  ES(]  ""*&)  '     Xd>^  (44) 

0.  otherwise 


Ea^l^-fc)    >     **>¥  (45) 

0,  otherwise 


where  adH  and  aaH  are  the  Hydrogen  radii  for  donors  and  acceptors,  respectively.  Equa- 
tions (44-^45)  are  the  analytical  solution  for  the  binding  energies  of  donor  and  acceptor 
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states  under  the  Hulthen  potential,  which  is  an  excellent  approximation  to  the  screened 
Coulomb  potential.  So,  the  Eqs.  (35-41)  are  solved  self-consistently  with  static  screening 
included  until  the  binding  energies  converge. 

Electronic  Mobility  due  to  Ionized  Impurity  Scatterings  [96] 

After  the  chemical  potential  is  determined  self-consistently,  the  screening  length 
under  the  generalized  Debye  theory  is  given  by  Eq.  (43).  The  mobility  due  to  ionized 
impurity  scatterings  is  then  computed  as  described  below,  including  the  static  screening 
effect  on  the  scattering  potential  of  the  ionized  impurities  from  other  charged  entities. 

In  order  to  obtain  the  electronic  mobility  due  to  ionized  impurity  scatterings,  the 
BKE  is  solved  under  the  relaxation  time  approximation.  As  generally  known,  rare 
analytical  solutions  are  found  for  BKE  [37].  It  is  worth  mentioning  that  when  dealing 
with  polar  LO  phonon  scatterings  at  low  temperature,  which  are  inelastic,  the  relaxation 
time  approximation  is  invalid.  A  variational  method  may  be  applied  to  compute  transport 
properties  [49].  For  materials  where  elastic  scatterings  dominates,  the  relaxation  time 
approximation  can  be  used  for  calculations  of  transport  parameters.  As  for  ionized 
impurity  scatterings,  the  scattering  events  are  elastic,  so  use  of  the  relaxation  time 
approximation  is  allowed.  In  order  to  compute  transport  coefficients,  it  is  needed  to 
find  the  relaxation  time,  which  is  directly  related  to  the  scattering  rates  for  different 
elastic  mechanisms. 

For  n-type  ZnSe,  only  electronic  scatterings  within  the  conduction  band  need  to  be 
considered,  which  means  both  minority  carrier  and  impurity  conductions  are  ignored.  The 
scattering  rate,  w(k,k'j,  of  an  electron  being  scattered  from  initial  state  \ck  >  to  final 
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state  \ck'  >  induced  by  interaction  Utnt(r)  can  be  obtained  using  the  Fermi's  golden  rule, 

w{lQ)^\H^8{Ei-E9), 

V        '        *  (46) 

Hkk>  =  Jd3rrck^rlUint(f)wJf), 

where  0  r{r)  is  the  Bloch  wavefunction  for  a  conduction  electron  in  k  state, 

1     1-  (47) 

=  4=e  «  r(f)  . 
V  is  the  crystal  volume.  Parabolicity  is  assumed  for  E%,  and  intervalley  scatterings  are 
not  considered.  This  is  reasonable  owing  to  large  energy  separation  and  weak  interaction 
at  large  wavevectors.  Also  neglected  is  the  periodic  part  of  the  Bloch  wavefunctions, 
which  is  equivalent  to  the  assumption  that  the  overlap  integral  between  the  periodic  parts 
of  the  initial  and  final  states  is  unity.  For  ionized  impurity  scatterings,  the  interaction 
Uint{r)  =  eV'(r)  ,  where  V(r)  is  given  as  in  Eq.  (81)  with  Z,e,  being  the  amount  of 
charge  on  an  ionized  impurity,  and  the  interaction  matrix  element  is  computed  to  be 

**        snr  (|)  +  a'2 
where  8  -  cos-1  (k  ■  k'j  is  the  angle  between  wavevectors  k  and  fc*,  and 

^=7T&<         a  =  2k\D.  (49) 

Then,  the  relaxation  time,  or  the  linear  momentum  relaxation  time,  which  is  favored 
by  large  angle  scatterings,  is  computed  as 

-  =  —r  fd3k'  W(k,  k')(l-  cos 6) .  (50) 

T       (2tt)37  V        /v 

Since  electrons  follow  FD  statistics  and  may  have  a  broad  energy  distribution,  the  above 

relaxation  time  normally  depends  upon  the  energy  of  the  scattered  electron.    A  proper 
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averaging  procedure1  is  therefore  needed  for  r  [97].    It  is  shown  that  the  appropriate 

procedure  is  to  average  over  r  weighted  with  the  energy  of  the  scattered  electron, 

<Et> 

<<  r  >>= ,  (5 

<  E  > 


where 


<Et>=    ldEg(E)fFD(E)Er, 
<E>=  JdEg(E)fFD(E)E,  (52) 


Ifd(E) 


exP  (itiT1)  + l 
where  g(E)  is  the  density  of  states  and  proportional  to  \/E  for  a  parabolic  band. 

For  screened  Coulomb  potential  of  ionized  impurities  of  type  i  with  density  nz,  an 

analytical  solution  can  be  found  for  -r,  assuming  independent  scatterings  from  different 

'i 

ionized  impurities,  and  the  result  is 

1  fZtele\2  1  In  (1+ a2) -if 


(53) 


which  is  exactly  the  same  as  derived  by  classical  treatment.  Adopting  the  £  £■  additive 
scheme  for  the  effective  p,  the  relaxation  time  is  found  to  be 


e2  Ei 


t  =  l6irJ2m*me— — -=■ —  (54) 

t 
where  the  summation  over  i  accounts  for  different  types  of  ionized  impurities  and  it  is 

assumed  that  all  Zx  equals  unity.  Finally  the  averaged  relaxation  time  is  found  by  using 

Eq.   (51)  and  the  associated  drift  mobility  is  given  by 

fid-=  ~—  «t»  .  (55) 

m:me 


1.  It  is  assumed  that  the  averaging  procedure  can  be  extended  to  quantum  systems  under  low 
electric  field. 
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The  difference  between  the  drift  and  Hall  mobility  is  studied  so  as  to  understand  the 
effect  of  statistics  on  the  two  mobilities.  For  the  known  E?  dependence  of  r,  the  two 
mobilities,  in  MB  distribution,  differ  by  a  constant  ratio  of  1.93  [97].  The  Hall  mobility 
is  derived,  assuming  that  electric  and  magnetic  fields  perturb  the  distribution  function 
independently  so  that  their  contributions  are  additive  to  the  unperturbed  distribution 
function,  as  follows 

<<t2>>  rsx\ 

HH  = Yp'd  '  ^    ' 

«  r  >> 

For  FD  distribution,  the  ratio  of  1.93  is  changed  because  of  the  statistical  degeneracy 
of  the  system.  In  most  cases  this  ratio  tends  to  be  close  to  unity,  which  is  verified  in 
this  research. 


CHAPTER  15 
RESULTS  AND  ANALYSIS 


Three  parametric  series  have  been  investigated:  (A)  fixing  the  donor  concentration 
Nd  and  increasing  the  acceptor  Na  till  0.5A^  to  study  the  effect  of  compensation;  (B) 
increasing  Nd  with  Na  set  to  zero;  (C)  increasing  both  concentrations  with  Nd  =  2.Ya. 
The  temperature  range  is  the  same  for  all  series  but  impurity  concentrations  vary  as 
above.  Series  A  is  studied  so  as  to  understand  the  transition  from  Series  B  to  Series 
C.  Series  C  is  for  compensated  systems  and  Series  B  is  for  uncompensated  systems. 
In  order  to  avoid  the  impurity  band  issue,  Nd  is  chosen  to  be  1017cm~3  for  Series  A, 
as  compared  to  the  critical  value  given  in  Eq.  (34').  The  experimental  variable  range 
is  qualitatively  covered  by  Series  C.  The  parametric  space  is  defined  in  this  study  as 
the  parameter  ranges  covered  by  the  three  series.  Theoretical  results  are  compared  with 
the  experiments  [87]  in  this  chapter  following  their  own  exposition.  Table  8  lists  the 
parameters  used  in  this  study. 

There  are  two  major  physical  processes  which  are  instrumental  in  understanding  all 
the  results  in  this  research:  (1)  compensation  by  acceptors  and  (2)  thermal  ionization1 
of  donors.  The  first  process  dominates  at  low  temperature  and  the  second  one  becomes 
important  and  dominant  when  the  donor  level  is  reduced  small  enough  to  generate  enough 
free  electrons  in  the  conduction  band  even  at  low  temperature.  Static  screening  is  en- 
hanced due  to  compensation  and  strongly  affects  the  second  process.  In  n-ZnSe,  acceptors 
play  no  direct  role  in  static  screening  except  compensation  and  scattering.    Because  it 


1.  Thermal  excitation  is  used  when  referring  to  electrons. 
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87 


is  assumed  that  acceptors  are  fully  ionized,  there  is  negligible  hole  concentration  in  the 
system  as  well.  Therefore,  donors  and  electrons  are  the  focuses  of  discussion  as  fellows 


Table  8:  Input  Parameters  Used  in  the  Calculations  of  Part  IV. 


Series  Parameter 

Values 

T(K) 

5 

8 

15 

17.5 

20 

22.5 

25 

28 

32  5 

38 

43 

50 

63.5 

70 

85 

100 

125 

150 

200 

250 

300 

A 

Nd  (cm-3) 

1017 

Na  (cm-3) 

10° 

101 

102 

103 

104 

105 

106 

107 

108 

109 

1010 

10u 

1012 

2.5  X 
10" 

5  x  10" 

10" 

2.5  X 
1013 

5x10" 

1014 

2.5  x 
1014 

5xl014 

10" 

2.5  x 
10" 

5  X  10" 

1016 

2.5  X 
10" 

5  X  1016 

B 

Nd  (cm-3) 

1013 

1014 

10" 

10" 

10" 

1018 

10" 

3  X  10" 

Na  (cm'3) 

0.0 

C 

Nd  (cm-3) 

10" 

1014 

10" 

10" 

10" 

10" 

10" 

3  X  10" 

Na  (cm-3) 



5  X  1012 

5  X  10" 

5X1014 

5x  10" 

5X  10" 

5x  1017 

5  X  10" 

1.5  X 

10" 

Series  A 


Series  A  is  designed  to  understand  the  interplay  of  both  major  processes  and  the 
interplay  of  screenings  by  ionized  donors  and  electrons.  It  is  interesting  to  find  that 
the  static  screening  effect  induces  a  first-order  phase  transition  with  temperature  for 
uncompensated  systems  when  the  donor  concentration  is  higher  than  1017cm""3  and  the 
transition  disappears  at  high  compensation. 
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Shown  in  Fig.  21  is  the  derivative  of  the  ionized  donor  concentration  with  respect 
to  the  chemical,  potential  versus  temperature.  The  derivative  quantitatively  indicates  the 
screening  capability  of  donors,  as  exhibited  in  Eq.  (43).  It  is  helpful  to  understand  the 
result  by  looking  at  the  uncompensated  system  at  low  temperature.  Firstly,  NJ  »  n 
holds  approximately.  Also,  it  is  found  from  Eq.  (92)  that  |-^-j  w  /3n  since  ionization 


r+ 


V  ' 

ratio  :jf-  is  low.  Furthermore,  the  electron  concentration  is  given  by 

n  =  NtefiE'  ,  (57) 

where  all  the  quantities  are  defined  in  Eq.  (38-41),  and  Ef,  the  Fermi  energy,  can  be 
assumed  to  be  -13meV  because  it  is  usually  pinned  at  the  middle  between  the  donor 
level  and  the  bottom  of  the  conduction  band.  So,  it  is  possible  to  estimate  the  Debye 
length  at  25  K  and  the  reduction  of  the  donor  level  as  a  first-order  approximation.  The 
Debye  length  is  about  300  A  and  the  reduction  of  the  donor  level  is  about  10  percent. 
Therefore,  static  screening  due  to  thermal  excitation  influences  the  system  critically,  given 
the  highly  nonlinear,  self-consistent  description  of  the  system.  The  transition  at  26.5  K 
is  a  near-metal-insulator  transition,  as  checked  with  the  following  results  of  the  chemical 
potential  (see  Fig.  26)  and  electron  concentration  (see  Fig.  27),  and  it  is  because  of  the 
static  screening  effect  due  to  thermal  excitation.  Discussion  of  the  nature  of  the  transition 
at  26.5  K  will  be  attempted  in  the  following  chapter.  For  a  compensated  system,  as 
acceptors  are  added  to  the  system,  electrons  in  donor  states  will  transfer  to  the  lower 
energy  states,  or  acceptor  states.  As  a  result,  the  acceptors  are  fully  ionized  presumably 
and  donor  states  are  partially  filled  so  that  they  will  participate  in  the  screening  piocess. 
At  zero  temperature,  iV+  =  iVa  so  that  |-^-|  =  ^Mi^zMt  according  to  Eq.  (92), 
is  finite,  which  means  finite  screening  capability.  Compensation  has  the  best  effect  on 
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screening  when  half  of  the  donors  are  ionized  according  to  Eq.  (92).  The  effect  of  the 
temperature  on  screening  is  twofold:  thermally  generating  electrons  in  the  conduction 
band  and  reducing  screening  through  the  (3  factor,  as  appeared  in  Eqs.  (88)  and  (92).  At 
low  temperature  and  low  compensation,  thermal  excitation  is  more  important  than  the  (3 
factor  effect,  as  analyzed  above.  When  temperature  rises,  the  static  screening  effect  is 
gradually  reduced  due  to  the  (3  factor. 

The  result  for  the  derivative  of  the  free  electron  concentration  with  respect  to  the 
chemical  potential,  as  shown  in  Fig.  22,  can  be  understood  along  the  same  line.  A  twist, 
however,  is  the  generation  of  electrons  in  the  conduction  band  at  iow  temperature.  Due 
to  the  fact  that,  the  donor  concentration,  A^,  is  high  enough,  1017cra~3,  a  MI  transition 
is  possible  for  a  highly  compensated  system.  By  MI  transition  it  is  meant  that  the 
chemical  potential  falls  into  the  conduction  band  since  the  bottom  of  the  conduction 
band  defines  the  mobility  edge  in  the  present  theory  When  a  MI  transition  occurs,  donor 
impurities  are  ionized  thermally  and  part  of  their  electrons  become  free  carriers  in  the 
conduction  band  Further  discussion  about  this  simple  picture  of  MI  transition  follows  in 
the  next  chapter.  So,  due  to  the  MI  transition  in  a  highly  compensated  system,  electrons 
appear  in  the  conduction  band  at  low  temperature.  And,  the  system  exhibits  the  behavior 
of  degenerate  electrons,  which  means  finite  screening  capability  at  low  temperature, 
according  to  \jj*-\  -  Bn  and  a  constant  n,  and  the  (3  factor  dominates.  However,  for 
systems  of  low  compensation  at  low  temperature,  free  electrons  arc  thermally  generated 
so  that,  n  follows  Eq.  (57)  and  the  screening  capability  grows  exponentially  before  the 
MI  phase  transition. 

It  is  straight-forward  to  deduce,  from  the  generalized  Debye  theory,  that  under  first- 
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order  approximation  fully  occupied  or  unoccupied  bands  does  not  screen  statically.  The 
Debye  length  of  a  rc-ZnSe  system  indicates,  inversely,  its  screening  capability.  As 
clearly  shown  in  Fig.  23,  when  the  system  has  no  acceptors,  the  Debye  length  goes 
to  infinity  at  low  temperature.  When  acceptors  are  added  to  the  uncompensated  n-ZnSe 
bulk  material,  according  to  the  present  theory,  all  of  them  compensate  donors  so  that 
donor  states  become  susceptible  to  perturbation  of  external  potential.  As  a  resuit,  the 
Debye  length  begins  to  decrease  and  reaches  minimum  at  5  x  1015cm-3  of  the  acceptor 
concentration.  Above  this  acceptor  concentration,  the  system  is  metallic  and  the  screening 
is  weakly  temperature-dependent  thanks  to  the  /?  factor  effect  This  critical  concentration 
is  defined  as  the  critical  acceptor  concentration  for  the  convenience  of  later  discussion. 
As  temperature  increases  for  uncompensated  systems,  the  Debye  length  decreases  and 
makes  a  transition  at  26.5  K  when  thermal  ionization  of  the  donor  impurities  critically 
enhances  the  screening  capability  of  donor  impurities  and  electrons.  This  transition 
temperature  is  defined  as  the  critical  temperature  for  later  discussion  although  it  is  not 
well-defined  at  high  acceptor  concentration.  Above  the  critical  temperature,  the  Debye 
length  increases  slightly  because  of  the  /3  factor  effect.  Over-doping  with  acceptors 
beyond  the  critical  acceptor  concentration  only  affect  the  screening  capability  of  the 
system  barely,  which  is  reasonable  because  of  the  maximum  in  the  screening  capability 
approximately  at  half  compensation. 

The  ionized  donor  concentration  is  important  in  determining  the  mobilities  of  free 
carriers  at  low  temperature  and  it  is  affected  by  compensation  and  thermal  excitation. 
In  the  present  theory,  all  acceptors  are  compensated,  or  ionized.  As  a  consequence, 
addition  of  acceptors  enhances  the  screening  capability  of  donors,  which  in  return  reduces 
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the  binding  energy  of  donors,  ^s  clearly  indicated  in  Fig.  24,  increase  of  acceptor 
concentration  reduces  the  binding  energy  of  donors  appreciably  from  their  initial  value 
E°d  -  27.3 meV,  as  given  in  Table  7.  The  result  in  Fig.  24  can  be  characterized  by  two 
values  as  before:  the  critical  acceptor  concentration  and  the  critical  temperature.  Below 
the  critical  temperature,  donor  level  is  reduced  with  acceptor  concentration  because  of  the 
compensation  effect.  The  (3  factor  effect  allows  the  donor  level  to  recover  appreciably 
due  to  weakening  in  the  screening  capability  of  the  system  Interestingly,  the  donor 
level  can  recover  somewhat  from  the  near-MI  transition  above  the  critical  temperature. 
More  surprisingly,  the  data  show  re-entrant  behavior  with  temperature  at  the  acceptor 
concentration  of  5  x  1015cm~3.  On  the  other  hand,  the  ionized  donor  concentration 
basically  follows  the  Debye  length  in  behavior.  Above  the  critical  acceptor  concentration, 
it  increases  slightly  with  temperature  due  to  thermal  ionization.  There  is  a  drastic  increase 
in  the  ionized  donor  concentration  in  the  transition  region. 

Due  to  the  reduction  of  donor  levels  by  screening  effect,  the  chemical  potential 
shows  a  sign  of  a  metal-insulator  phase  transition  at  a  critical  temperature  of  26.5  K  for 
uncompensated  systems  (25/  =  -0.094  me V  at  26.5  K),  but  retreats  gradually  because 
of  lack  of  screening  capability  to  sustain  the  transition  at  the  donor  concentration  of 
1017cm-3  for  Series  A.  Above  the  critical  temperature,  the  chemical  potential  begins  to 
exhibit  intrinsic  behavior  of  the  bulk  material  and  decreases  towards  the  middle  of  the 
bandgap.  Thermodynamically,  it  is  easy  to  see  that  below  the  critical  temperature  donor 
level  is  the  dominant  energy  scale  but  static  screening  effect  strongly  modifies  it  so  as 
to  cause  the  system  to  show  sign  of  phase  transition.  Above  the  critical  temperature, 
the  bandgap  of  the  host  material  kicks  in  as  the  smaller  energy  scale,  i.e.,  the  donor 
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level,  is  submerged  by  thermal  fluctuations.  Along  the  same  line  as  before,  the  critical 
acceptor  concentration  distinctively  classifies  the  result  in  Fig.  26  into  two  regions:  above 
the  critical  concentration,  the  system  is  metallic  (Ef  =  0.6  meV  at  6  K  at  the  critical 
concentration)  and  thermal  excitation  drives  the  chemical  potential  towards  the  middle 
of  the  bandgap;  below  the  critical  concentration,  the  chemical  potential  is  pinned  around 
half  of  the  donor  level  at  low  temperature  and  then  the  chemical  potential  experiences  a 
near-MI  transition  because  of  the  decrease  of  the  donor  energy  due  to  screening  effect. 

Free  electrons  in  the  conduction  band  are  generated  by  thermally  ionizing  donors. 
So,  the  free  electron  population  reflects  the  thermodynamical  activity  of  donor  impurities. 
At  low  acceptor  concentration,  as  shown  in  Fig.  27,  there  is  an  abrupt  change  in  the 
free  electron  density  at  26.5  K  and  the  electron  density  is  altered  by  about  three  orders 
of  magnitude.  This  is  a  clear  manifestation  of  the  near-MI  transition  due  to  the  static 
screening  effect.  Below  the  critical  temperature,  the  density  follows  the  behavior  of 
thermal  excitation  over  a  single  energy  gap,  which  shows  a  linear  dependence  in  the 
logn-jf  plot.  This  behavior  is  known  as  the  freeze-out  behavior.  Saturation  is  nearly 
reached  above  the  critical  temperature  where  the  electron  density  is  almost  independent  of 
temperature.  As  acceptor  density  increases,  the  critical  temperature  is  reduced  gradually 
until  a  MI  transition  is  reached.  Around  this  transition,  re-entrant  behavior  has  been  found, 
which  presumably  has  to  do  with  the  highly  nonlinearity  near  the  transition  region.  In 
the  metallic  regime,  the  electron  density  increases  gradually  with  temperature  though  it 
can  become  twice  as  much  from  77  K  to  300  K. 

Drift  mobility  and  Hall  mobility  limits  how  fast  a  semiconductor  can  operate.  When 
semiconductors  are  doped,  it  is  desired  to  have  as  high  a  mobility  as  possible.   At  low 
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temperature,  carrier  scatterings  with  charged  impurities  essentially  determine  the  mobil- 
ities, though  scatterings  with  phonons  dominate  at  room  temperature.  Therefore,  result 
for  mobilities  below  50  K  is  reported.  Figure  28  shows  the  Hall  mobility  as  a  function  of 
temperature  and  acceptor  concentration  and  the  drift  mobility  is  characteristically  similar. 
Below  the  critical  acceptor  concentration  and  the  critical  temperature,  the  Hall  mobility 
drops  about  ten  orders  of  magnitude,  but  changes  barely  in  other  regions.   This  result 
is  in  sharp  contrast  with  predictions  by  the  well-known  CW  [89,  96]  and  BH  [88,  96] 
theories.   Both  theories  predict  no  drastic  change  in  the  mobilities.   The  difference  re- 
sults, obviously,  from  whether  the  static  screening  effect  and  the  chemical  potential  are 
treated  self-consistently.  According  to  Eq.  (54-55),  drift  mobility  depends  upon  ionized 
impurity  concentrations  inversely  and  screening  weakly.  Thus,  in  regions  where  ionized 
impurity  concentrations  vary  rapidly,  the  mobility  is  controlled  mainly  by  ionized  im- 
purity concentrations,  which  is  supported  by  Fig.    25  in  the  low  temperature  and  low 
acceptor  concentration  region.   In  other  regions  where  ionized  impurity  concentrations 
vary  inappreciably,  screening  effect  will  play  a  major  role,  though  it  does  not  affect  the 
mobilities  strongly. 

As  a  note,  conductivity  or  resistivity  of  the  doped  material,  in  case  of  doping  a 
semiconductor  to  make  ohmic  contact,  is  the  critical  parameter  and  the  conductivity  by 
electrons  depends  upon  the  product  of  the  electron  concentration  and  the  drift  mobility. 

Series  B  and  C 

As  clearly  exhibited  in  Figs.  29-37,  highly  compensated  and  uncompensated  semi- 
conductors show  characteristically  different  behaviors  owing  to  enhanced  screening  effect 
in  the  compensated  materials  (Series  C)  by  compensation.    As  suggested  in  Eq.    (77), 
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any  deviation  from  full  or  zero  occupation  of  an  energy  state  will  enhance  its  screening 
capability  until  half  occupancy  For  n-type  semiconductors,  the  addition  of  acceptors 
is  instrumental  to  this  screening-capability-enhancing  process  for  donors.  This  process 
makes  a  fundamental  difference  in  the  two  series  at  low  temperature.  However,  as  tem- 
perature goes  up,  thermal  ionization  of  donor  states  will  smear  out  this  difference,  as 
already  exhibited  in  the  results  for  Series  A  and  shown  in  the  following  results. 

At  low  temperature,  donor  ionization  is  determined  by  compensation  in  Series  C  but 
by  thermal  ionization  in  Series  B.  The  ionized  donor  concentration  equals  the  accep- 
tor concentration  at  zero  temperature  for  Series  C  when  the  donor  concentration  is  low. 
However,  when  the  donor  concentration  is  over  101' cm-3,  the  system  becomes  metal- 
lic and  the  ionized  donor  concentration  saturates  at  all  temperature.  The  ionized  donor 
concentration  does  not  change  much  in  magnitude  in  Series  C  since  half  of  the  donors 
have  already  been  ionized,  due  to  compensation,  even  at  low  donor  concentration.  As. 
temperature  increases,  the  transitional  behavior  in  donor  ionization  in  Series  B  sets  in 
thanks  to  thermal  ionization  and  the  critical  temperature  decreases  with  donor  concen- 
tration because  screening  capability  gets  stronger  since  the  capability  is  proportional  to 
the  donor  concentration,  as  exhibited  in  Fig.  30.  The  MI  transition  with  temperature  in 
Series  B  is  well-defined  above  1017cm~3  of  donor  concentration,  but  in  Series  C  no  such 
abrupt  transitions  are  exhibited  with  temperature  since  compensation  dominates  similar 
to  Series  A  in  the  previous  section.  Above  the  critical  temperature,  as  shown  in  Fig. 
29,  ionized  donor  concentration  is  almost  the  same  for  the  two  series  and  increases  with 
donor  concentration.  For  both  series,  donor  ionization  by  thermal  excitation  tends  to 
saturate  towards  high  temperature  for  non-metallic  systems,  as  expected. 
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As  shown  in  Fig.  30,  at  low  donor  concentration,  for  both  series,  the  screening  effect 

is  weak  as  the  Debye  length  is  over  1,000  A.    As  donor  concentration  increases,  the 

_  i_ 
screening  capability  increases  roughly  according  to  Nd  2  for  Series  C,  and  this  behavior 

holds  for  Series  B  above  the  critical  temperature.  This  behavior  can  be  seen  from  Eq.  89 

and  the  fact  that  n  ex  Nd  when  Nf  varies  slowly.  In  all  temperature  ranges,  for  Series 

C.  the  temperature  dependence  is  roughly  v'T,  as  predicted  by  the  (3  factor  effect.  For 

Series  B,  the  temperature  dependence  features  a  critical  behavior  at  the  transition  and 

weak  dependence  above  it  at  high  donor  concentrations. 

As  shown  in  Fig.  31,  the  dependence  of  the  Debye  length  on  carrier  concentration  for 
Series  B  does  go  roughly  as  the  inverse-square-root  prediction  for  classical  electron  gas 
below  i018cm-3  of  the  election  density.  The  reason  is  because  that  ionized  donors 
contribute  to  screening,  according  to  Figs.  21  and  22,  with  the  same  dependence. 
Therefore,  the  dependence  follows  the  prediction  for  classical  electron  gas  at  low  electron 
density.  The  deviation  at  high  electron  density  in  Series  B  is  due  to  the  0  factor  effect 
above  the  transition  region.  The  Debye  length  for  Series  C  is  weakly  dependent  on 
the  electron  concentration  for  systems  with  low  donor  concentration  and  the  reason  is 
because  compensation  dominates  over  thermal  excitation  of  electrons  in  contributing  to 
the  screening  capability.  The  deviation  at  high  electron  density  in  this  series  is  due  to 
the  same  reason  as  in  Series  B.  However,  positive  dependence  of  a  power  of  .33  at  high 
concentration  is  shown  and  explanation  is  the  /?  factor  effect  because  the  Debye  length 
is  reduced  by  the  effect  while  the  electron  concentration  does  not  change  noticeably  in 
these  metallic  systems,  as  shown  in  Fig.    34. 

Donor  level  directly  relates  to  the  screening  capability  of  the  system  according  to 
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Eq.    (44-45),  as  reiterated  below, 

Ei  =  lE}(l-{&f,   ad>4  (58) 


V' 


otherwise 


As  estimated  in  Table  7,  the  hydrogenic  radius,  adH,  of  a  donor  in  ZnSe  is  28.6  A. 
The  present  theory  assumes  that  Er(  is  reduced  to  zero  if  the  Debye  length  is  shorter  than 
half  of  the  donor  hydrogenic  radius.  As  shown  in  Fig.  32,  donor  level  shrinks  in  Series 
B  at  low  temperature  and  reaches  2ero  when  donor  concentration  exceeds  101  cm. 
Around  this  donor  concentration,  the  donor  level  tends  to  recover  as  the  temperature 
increases,  which  is  because  the  /?  factor  effect  decreases  the  screening  capability  of 
the  system.  The  donor  level  for  Series  C  increases  smoothly  with  temperature  and 
decreases  with  donor  concentration.  It  follows  the  behavior  of  the  Debye  length:  (1)  the 
screening  increases  with  donor  concentration;  (2)  at  low  donor  concentration  the  thermal 
excitation  contributes  more  appreciably  to  screening  than  at  high  donor  concentration, 
due  to  compensation;  (3)  the  donor  level  shrinks  to  zero  above  1018cm-3  of  donor 
concentration;  (4)  the  /3  factor  effect  reduces  screening  with  temperature. 

The  chemical  potential  defines  the  state  of  the  system:  at  low  temperature,  it  becomes 
positive  at  MI  transition;  at  high  temperature,  because  the  intrinsic  behavior  sets  in,  it 
reduces  towards  the  middle  of  the  bandgap  even  though  the  system  is  still  metallic,  as 
indicated  by  the  free  electron  density!  The  chemical  potential  behaves  quite  differently  for 
Series  B  and  C  as  shown  in  Fig.  33  The  main  difference  exists  in  the  donor  concentration 
dependence,  which  is  a  reflection  of  the  acceptor  concentration  in  the  series,  as  given 
in  Table  8.  Series  B  displays  MI  transitional  behavior  above  donor  concentration  of 
1017cm-3  with  temperature,  induced  by  the  screening  effect  due  to  thermal  excitation, 
and  the  critical  temperature  decreases  with  donor  concentration  because  of  increase  in 
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the  screening  capability  of  the  system.  On  the  other  hand,  the  chemical  potential  in  Series 
C  gradually  increases  and  finally  merges  into  the  conduction  band  at  low  temperature. 
This  leads  to  a  MI  transition.  The  critical  donor  concentration  is  about  10  cm"3, 
roughly  comparable  with  the  critical  donor  concentration  as  predicted  in  Eq.  (34)  tor 
the  Anderson  Ml  transition.  Tiie  implication  of  this  coincidence  will  be  discussed  in  the 
following  chapter.  In  comparison,  Series  C  has  a  weaker  temperature  dependence  than 
Series  B  because  of  compensation  (recall  that  the  acceptor  concentration  is  zero  in  Series 
B  and  half  of  the  donor  concentration  in  Series  C).  At  high  temperature,  the  intrinsic 
behavior  of  the  bulk  semiconductor  dominates,  as  the  valence  bands  pull  the  chemical 
potential  down  Generally,  the  chemical  potential  is  lower  in  Series  C  than  in  Series  B, 
as  acceptors  pull  the  chemical  potential  down. 

The  electron  concentration  measures  the  doping  efficiency  of  a  doped  semiconductor 
and  affects  the  conductivity  of  the  material.  For  donor-only  semiconductors  of  low  donor 
concentration,  the  well-understood  behavior  of  the  electron  concentration  features  three 
regimes  in  temperature:  (1)  thermal  excitation  over  donor  level  at  low  temperature,  (2) 
saturation  of  electron  concentration  at  temperature  comparable  to  the  donor  level  and 
(3)  thermal  excitation  over  bandgap  at  temperature  comparable  to  the  bandgap.  The 
result  for  Series  B  exhibits  this  typical  behavior  at  low  donor  concentration.  Obviously, 
room  temperature,  25.9  meV,  is  far  too  low,  compared  to  2,7  eV,  the  bandgap  of  ZnSe, 
to  see  the  high  temperature  regime.  The  slope  in  the  logn-4  plot  for  the  freeze-out 
regime  is  measured  at  -0.1462  K,  which  yield  a  value  of  25  2  meV  for  E&  according 
to  a  ex  exp  ( ~  ajftr )  ■  The  MI  transition  with  temperature  takes  place  above  a  donor 
concentration  of  1017 cm"3.  On  the  other  hand,  Series  C  shows  a  transitional  behavior 
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at  a  donor  concentration  of  1017 cm'3.  In  addition.  Series  B  always  generates  an  order 
of  magnitude  more  free  carriers  with  the  same  donor  concentration  than  Series  C  in  the 
metallic  region.  Therefore,  the  doping  efficiency  for  Series  B  is  higher  than  for  Series  C. 

The  Hall  mobility  is  another  important  property,  other  than  free  electron  concen- 
tration, that  affects  the  conductivity  of  a  doped  semiconductor.  Figure  35  shows  that 
Hall  mobility  for  Series  B  and  C  differs  drastically  below  the  critical  temperature.  The 
ionized  impurity  population  increases  from  zero  in  Series  B  and  ionized  donors  and  ac- 
ceptors exist  in  Series  C  even  at  zero  temperature,  thanks  to  compensation.  In  Series 
B,  below  the  critical  temperature,  mobility  drops  exponentially  with  temperature,  due  to 
exponential  increase  in  ionized  donor  concentration.  This  decrease  slows  down  above  the 
critical  temperature  in  the  metallic  phase  because  of  the  strong  screening.  In  Series  C,  the 
mobility  increases  at  low  donor  concentration  but  decreases  at  high  donor  concentration 
A  reversal  occurs  at  1016cm~3.  Equations  (54-55)  tell  that  mobility  depends  inversely 
upon  ionized  impurity  concentration  and  weakly  upon  the  screening,  if  the  screening 
effect  behaves  smoothly,  as  in  Series  C.  In  order  to  analyze  the  temperature  dependence 
of  the  mobility,  it  is  necessary  to  look  at  Eq.  (54-55)  carefully.  The  screening  effect  is 
weak  when  donor  concentration  is  low.  Owing  to  compensation,  the  ionized  impurity 
concentration  does  not  change  much,  as  shown  in  Fig.  29,  so  that  screening  effect  di- 
rectly influences  mobility.  At  low  temperature  the  relaxation  time,  according  to  Eq.  (54), 
is  proportional  to  E*.  Since  the  screening  effect  is  weak,  the  mobility  increases  with 
temperature,  as  the  classical  electron  gas  has  <  E  >-  \ksT.  On  the  contrary,  at  high 
donor  concentration  and  low  temperature,  the  screening  effect  is  strong  and  the.  a  factor 

d 

in  Eq.  (54)  gets  small  (a  =  0.33  if  A/j  =  -f-  and  k  is  estimated  from  classical  electron 
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_    i 
gas  with  T  =  25  K)  so  the  energy  dependence  of  the    relaxation  time  changes  to  E   3 

(an  expansion  in  a  produces  an  E2  term  in  the  denominator  in  Eq.  (54).),  which  means 

that  mobility  decreases  with  temperature  instead. 

The  Hall  mobility  is  not  an  explicit,  function  of  the  electron  density  in  the  present 
theory  (since  carrier-carrier  scattering  is  not  included)  and  is  shown  in  Fig.  36.  Two 
regimes  can  be  classified  for  both  series:  (1)  the  low  electron  density  regime  and  (2)  the 
high  electron  density  regime.  At  low  electron  density.  Series  B  can  be  characterized  by 
a  well- defined  power  law  behavior  with  a  power  of  -.89,  and  above  the  MI  transition, 
a  different  power  of  approximate  -1.6  can  be  used  to  describe  the  relation.  Series  C 
varies  insensitively  with  carrier  concentration  at  low  electron  density  but  converges  to  a 
power  law  behavior  after  the  MI  transition. 

The  ratio  of  the  Hall  mobility  to  the  drift  mobility  indicates  the  statistics  of  free 
carriers.  For  a  MB  distribution,  various  theories  [96]  predict  approximately  1.9  for  the 
ratio,  and  for  FD  distribution,  it  is  predicted  to  be  about  unity.  As  shown  in  Fig.  37, 
the  ratio  is  bounded  between  roughly  1.0  and  1.9.  It  basically  follows  the  variation  of 
the  donor  level  shown  in  Fig.  32.  At  low  donor  concentration,  the  two  series  converges 
to  the  MB  value,  and  to  the  FD  value  at  high  donor  concentration  except  for  Series  B 
at  low  temperature.  At  low  temperature,  the  ratio  approaches  the  classical  limit  of  1.9 
in  Series  B,  but  there  is  no  convergence  in  Series  C  and  the  ratio  depends  upon  details 
of  the  static  screening  effect.  A  donor  concentration  of  101Tcm~3  is  the  critical  value 
where  systems  begin  to  fall  into  different  statistical  limit.  For  n-ZnSe,  the  degree  of 
degeneracy  of  a  system  must  be  considered,  following  Ruda  [49].    A  critical  electron 
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density  can  be  defined  as 
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n„  =  NcF±{0) 


(59) 


from  Eq.  (38),  which  indicates  that  FD  statistics  must  be  applied  to  a  system  with 
an  electron  density  larger  than  this  critical  value.  Thus,  ncr(300A')  =  1.34  x  10  cm"  , 
ncr(77A)  =  1.75  x  1017cm"3  and  ncr(25A)  =  3.24  x  1016cm-3.  This  last  value  crudely 
corresponds  to  the  critical  donor  concentration  found,  though  compensation  makes  the 
two  critical  values  differ. 

Experimental  Results 

The  experimental  data  by  Park  and  Simmons  [87]  are  listed  below  for  a  qualitative 
comparison  to  the  present  theory.  The  experiments  are  on  heavily  doped  and  interme- 
diately doped  n-ZnSe  semiconductors.  Table  9  enumerates  the  impurity  concentrations 
for  all  six  samples.  For  convenience,  the  ratio  of  acceptor  concentration  to  donor  con- 
centration is  also  calculated  in  the  table. 

The  chemical  potential  for  Sample  2-5  are  measured  and  Sample  5  is  metallic  as  it 
stays  within  the  conduction  band,  as  shown  in  Fig.   38,  and  its  free-electron  density  is 


Table  9:  Donor  and  Acceptor  Concentrations  in  the  Experiments  [87]. 


Sample 

Nd  (I016cm"3) 

Na  (I016cm-3) 

& 

1 

2.6 

0.9 

0.35 

2 

4.6 

1.4 

0.30 

3 

20.5 

7.7 

0.37 

4 

61.0 

35.0 

0.57 

5 

153.0 

103.0 

0.67 

6 

2370.0 

1510.0 

0.64 
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a  constant  and  high,  as  shown  in  Fig.  39.  Sample  4  shows  a  crossover  behavior  while 
the  other  two  samples  remains  semiconducting. 

The  free  electron  density  behaves  distinctively  in  different  phase:  in  the  metallic 
phase,  the  density  does  not  change  with  temperature,  however,  in  the  semiconducting 
phase,  freeze -out  and  saturation  regimes  are  distinguishable.  All  samples  have  been 
measured  for  their  free  electron  density,  and  Sample  1  and  2  display  semiconducting 
behavior,  Sample  5  and  6  exhibit  metallic  trait,  with  Sample  3  and  4  showing  deviation 
towards  metallic  characteristics  though  remaining  semiconducting,  as  shown  in  Fig.  39. 

Also  measured  for  all  the  samples  is  the  dependence  of  the  Hall  mobility  on  tempera- 
ture. Sample  5  and  6  show  metallic  behavior  in  a  heavily  doped  system,  i.e.,  independent 
of  temperature.  The  other  samples  show  temperature  dependence  similar  to  intermediately 
doped  systems.  At  high  temperature,  the  phonon  contribution  dominates  and  the  mobility 
decreases  with  temperature  as  the  magnitude  of  the  characteristic  wavevector  of  phonons 
increases.  However,  at  low  temperature,  the  phonon  population  drops  and  ionized  impuri- 
ties dominate  the  momentum-loss  processes  of  free  carriers.  The  momentum-loss  process 
due  to  ionized  impurity  scatterings  depends  upon  the  characteristic  energy  and  momen- 
tum of  the  carrier,  as  shown  in  Eq.  (54),  which  increases  with  temperature.  Hence,  the 
mobility  increases  with  temperature.  The  crossover  temperature  from  phonon  domination 
to  ionized  impurities  domination  increases  with  impurity  concentration. 

Qualitative  comparison  of  these  experimental  data  with  the  results  of  this  research 
follows.  Comparing  Figs.  38,  33,  39,  and  34,  similarities  are  seen  between  the 
experimental  data  and  the  results  for  Series  C  with  the  present  theory.  For  instance, 
between  Sample  4  and  the  calculations  with  donor  concentration  of  1017cm~"  and 
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1018cm~3  for  Series  C,  all  of  them  show  crossover  behavior  in  chemical  potential  and 
the  free  electron  density  varies  with  temperature  in  the  same  style.  However,  from 
the  limited  experimental  data  on  low  temperature  Hall  mobility,  the  disagreement  is 
overwhelming  between  experiments  and  theory  at  high  donor  concentration.  Not  only 
the  theoretical  prediction  is  many  orders  of  magnitude  larger  than  the  experimental  data, 
but  also  they  completely  differ  in  their  temperature  dependence.  The  possible  causes  for 
disagreements  will  be  discussed  in  the  next  chapter. 
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Figure  38:  Measured  chemical  potential  for  n-ZnSe  [87]. 
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Figure  39:  Measured  free-electron  density  for  n-ZnSe  [87]. 
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CHAPTER  16 
DISCUSSION 


First  of  all,  the  static  screening  effect  is  not  negligible1  in  semiconductors  even 
at.  a  donor  concentration  as  low  as  1015cm~3,  as  suggested  in  Fig.  30.  The  critical 
parameter  is  the  ratio  of  the  Debye  length  to  the  donor  hydrogenic  radius  for  the  n-type 
semiconductor.  For  III-V  and  II- VI  semiconductors,  the  donor  hydrogenic  radius,  as 
estimated  by  assuming  that  the  relative  effective  mass  is  0.1  and  the  relative  dielectric 
constant  is  70,  is  one  hundred  times  of  the  Bohr  radius,  according  to  adH  =  .529/c/m*. 
Thus,  a  Debye  length  with  a  value  of  thousands  A  will  modify  donor  level  appreciably, 
according  to  Eq.  (44-45).  This  means  that  impurities  can  change  the  physical  properties 
of  the  host  material  at  a  much  lower  concentration  than  the  critical  concentration  for 
metal-insulator  transition. 

The  static  screening  effect  modifies  material  parameters  and  may  induce  a  phase 
transition.  The  Debye  length  affects  the  impurity  levels.  Thermodynamic  equilibrium 
of  the  independent  particle  system,  as  considered  in  the  present  theory,  is  related  to 
compensation  and  thermal  excitations  which  are  directly  influenced  by  impurity  levels. 
The  two  physical  processes  determine  the  Debye  length.  This  nonlinear  enhancement  of 
screening  may  iead  to  critical  behavior,  as  shown  in  the  previous  chapter.  As  for  those 
controlling  factors  for  the  critical  temperature,  which  no  analytical  expression  has  been 
found,  calculations  have  been  done  concerning  the  functional  form  of  the  donor  level 
reduction,  i.e.,  Eq.   (44-45),  dielectric  constant,  and  the  initial  donor  level.  It  is  found 


1.  It  has  been  observed  that  the  donor  binding  energy  vanishes  for  impurity  concentrations  as 
low  as  1014cm-3  in  InSb  (see  reference  98). 
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that  change  of  the  functional  form  from  quadratic  to  linear  form  diminishes  the  transition, 
as  shown  in  Fig.  41.  The  dielectric  constant  k  is  artificially  varied  without  adjusting  E°[ 
accordingly  and  the  critical  temperature  diverges  when  lowering  n  and  saturates  towards 
about  22  K  when  increasing  k.  When  oniy  varying  E%;  since  £°  has  no  effect  on 
the  critical  temperature  as  expected  and  checked,  it  shows  that  the  temperature  roughly 
vanes  linearly.  Thus,  we  guess  that  the  critical  temperature  is  linearly  proportional  to 
the  donor  level  E\,  as  intuitively  expected. 

The  chemical  potential  varies  in  the  parametric  space  and  strongly  affects  the  screen- 
ing capability  of  the  doped  semiconductor.  The  screening  effect  depends,  primarily,  upon 
the  accessibility  of  those  charged  and  thermodynamically  active  particles,  which  is  deter- 
mined by  the  chemical  potential.  This  nonlinear  relationship  means  that  it  is  not  plausible 
to  assume  a.  constant  chemical  potential  when  dealing  with  charged  particles.  Therefore, 
self-consistency  is  enforced  in  studying  the  whole  parametric  space. 

Also,  the  chemical  potential  is  directly  related  to  the  generation  of  free  carriers,  and 
ihe  current-carrying  capability  of  the  doped  material.  Thus,  it  is  desired  to  increase  the 
efficiency  of  generating  free  charge  carriers  with  doping.  However,  adding  acceptors  to 
an  n-typc  semiconductor  results  in  both  increase  of  ionized  impurities  and  lowering  of 
the  chemical  potential.  As  a  consequence,  this  leads  to  decrease  in  mobility  and  free 
carrier  concentration.  Therefore,  it  is  impossible  to  dope  an  intrinsically  compensated 
semiconductor  efficiently,  unless  making  the  material  uncompensated. 

The  present  theory  shows,  partially  and  qualitatively,  agreement  behavior  with 
the  experimental  data.  Several  elements  hinder  attempt  of  quantitative  analysis  of 
the  experimental  data.    Firstly,  the  present  theory  does  not  have  any  fitting  material 
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parameters.  For  example,  the  donor  and  acceptor  levels  are  computed  self-consistently, 
not  as  fitting  parameters.  Furthermore.,  The  experiments  [87]  for  n-ZnSe,  which  is 
an  intrinsically  compensated  semiconductor,  do  not  have  reliable  data  for  impurity 
concentrations.  Last  but  not  the  least,  no  low-temperature  mobility  data  for  ZnSe  are 
available  for  comparison  where  the  ionized  impurity  contribution  predominates. 

Several  issues  [96]  may  affect  the  outcome  of  the  present  theory.  Discussed  here  are 
several  important  ones.  First  of  all,  the  generalized  Debye  theory  treats  static  screening 
perturbatively  which  definitely  overestimates  the  effect,  especially  at  high  donor  concen- 
trations and  MI  transition  region.  For  instance,  the  theoretical  mobilities  are  not  only 
three  orders  of  magnitude  higher  than  the  measured  values  but  also  show  different  temper- 
ature dependence  at  high  donor  concentration  (see  Figs.  35  and  40).  New  experimental 
data  (see  reference  28,  page  88)  on  high-quality,  CI  doped,  uncompensated  n-ZnSe/GaAs 
show  a  high  free-electron  density  up  to  1019cm~3,  semiconducting  behavior,  and  no  MI 
transition,  which  indicate  (1)  the  Debye  theory  does  overestimate  static  screening  and  (2) 
static  screening  is  important,  as  deduced  from  the  high  carrier  concentration.  But  this 
does  not  rule  out  possible  MI  transitions  in  other  materials  due  to  static  screening  effect 
[98].  For  static  electronic  screening  at  high  density,  many-body  theory,  such  as  random 
phase  approximation  (RPA),  may  be  adopted  [99].  Secondly,  in  the  present  theory,  ac- 
ceptors are  completely  compensated  and  no  correlations  at  high  concentrations  between 
impurities  are  considered.  Furthermore,  impurity  states  are  treated  by  the  hydrogenic 
model.  Only  ground  states  for  the  hydrogen-like  impurities  are  taken  into  account  in 
the  present  theory  but  excited  states  exist  in  real  systems  In  addition,  impurities  may 
be  doubly  occupied.  All  these  would  effectively  reduce  the  static  screening  due  to  im- 
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purities  at  low  temperature  since  the  impurities  may  not  be  ionized  when  perturbed,  or 
excited.  Another  assumption  made  in  this  reseauh  is  homogeneity  of  impurities  in  the 
system.   In  real  semiconductors,  inhomogeneity  exists,  which  will  appreciably  broaden 
DOS  distributions  for  both  electrons  and  impurities  at  high  concentrations.  This  would 
lead  to  increase  of  static  screening  as  the  mini-gap  between  the  donor  levels  (or  band)  and 
conduction  band  shrinks.  The  most  important  issue,  however,  is  that,  at  high  donor  con- 
centration, the  system  can  not  be  treated  as  an  ensemble  of  independent  impurities  since 
correlation  between  impurity  states  must  be  considered.  In  this  case,  many-body  theory 
is  needed  to  correctly  treat  the  system.  All  of  these  issues  may  influence  the  nature  of 
physics  involved  and/or  alter  the  results  of  this  research.  It  is  not  true,  however,  to  imply 
that  this  research  can  not  used  at  high  donor  concentration.  As  the  new  ZnSe/GaAs  data 
(see  reference  28,  page  88)  suggest,  static  screening  will  effectively  reduce  correlations 
in  the  system  at  high  donor  concentration  such  that  the  system  still  behaves  as  a  non- 
degenerate  semiconductor.  However,  considerations  do  have  to  be  taken  to  account  for 
the  overestimate  of  the  static  screening  effect,  especially  for  the  static-screening-induced 
MI  transition.  The  present  theory  does  agree  with  most  of  the  experimental  data.  Other 
than  these  major  issues,  the  following  considerations  may  .improve  the  present  theory:  (1) 
Born  approximation  may  be  replaced  by  more  accurate  partial  wave  phase-shift  analysis; 
(2)  bandstructure  effect  may  be  included  at  high  carrier  density;  and  (3)  multi-scattering 
events  may  be  considered  in  BKE. 

Finally,  this  research  is  restricted  to  ionized  impurity  scatterings.  However,  at  low 
temperature,  neutral  impurity  scatterings  may  become  important,  which  would  smear  out 
the  critical  behavior  in  mobility  in  certain  parametric  region.  Also,  as  temperature  rises, 
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phonon  contributions  compete  with  impurity  contributions  and  could  predominate  beyond 
50  K  for  semiconducting  systems.  Therefore,  full  consideration  of  all  major  scattering 
mechanisms  is  required  to  model  the  mobility  data  in  the  whole  temperature  range. 


CHAPTER  17 
CONCLUSIONS 


As  shown  in  this  research,  static  screening  effect  plays  an  important  role  in  determin- 
ing the  physical  properties  of  intermediately  to  heavily  doped  semiconductors,  especially 
compensated  ones.  For  n-ZnSe  is  an  intrinsically  compensated  semiconductor,  it  is  im- 
perative to  treat  the  static  screening  in  a  self-consistent  manner.  A  generalized  Debye 
theory  i S  adopted  for  this  purpose  and  the  analytical  result  of  the  Hulthen  potential  for 
impurity  levels  is  used.  The  drift  and  Hall  mobility  are  obtained  by  solving  the  Boltz- 
mann  kinetic  equation  under  the  relaxation  time  approximation  and  the  relaxation  time 
is  derived  quantum  mechanically  from  the  Born  approximation. 

Siatic  screening  effect  is  determined  by  compensation  and  thermal  ionization  of 
donor  imparities  in  «-ZnSe.  In  compensated  n-ZnSe,  at  low  temperature,  compensation 
creates  ionized  donors  and  acceptors,  which  effectively  enhances  (he  screening  capability 
and  thus  changes  the  characteristics  of  the  system.  And  as  temperature  rises,  thermal 
ionization,  influences  the  system  to  certain  degree  that  is  dependent  on  compensation 
ratio.  In  contrast,  in  uncompensated  n-ZnSe,  above  donor  concentration  of  approximate 
101TC77?"3,  a  well-defined  MI  transition  occurs  at  26.5  K,  but  is  smeared  out  when 
compensation  increases.  However,  no  such  transition  is  observed  experimentally,  which 
leads  the  author  to  conclude  that  the  Debye  theory  over-estimates  the  static  screening  at 
high  donor  concentration  when  the  screening  is  strong.  And  the  theory  is  applicable 
to  a  weakly  screened  system  with  its  Debye  length  at  least  five  times  larger  than  the 
donor  hydrogenic  radius.     Better  nonlinear  screening  theory  is  needed  at  high  donor 
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concentration.    Nevertheless,  the  results  in  this  research  are  completely  acceptable  at 
temperatures  and  donor  concentrations  away  from  the  transitional  region. 

It  is  found  that  static  screening  effect  greatly  affects  carrier  generation  and  mobility 
of  the  system  even  when  donor  concentration  is  far  below  the  critical  concentration  for 
metal-insulator  phase  transition.  Carrier  density  is  reduced  by  an  order  of  magnitude 
when  compensation  increases  at  high  donor  concentration  Therefore,  it  is  essential  to 
avoid  compensation  in  semiconductors  in  order  to  achieve  high  doping  efficiency  at  high 
donor  concentration.  Mobility  is  reduced  by  doping,  and  its  temperature  dependence,  in 
highly  compensated  systems,  changes  from  positive  to  negative  relationship. 

Finally,  this  research  partially  explains  the  experimental  data  by  Park  and  Simmons. 
The  chemical  potential  and  free  electron  density  data  are  well  understood  by  the  present 
theory  at  intermediate  doping.  Characteristically,  the  chemical  potential  decreases  with 
temperature,  and  increases  with  donor  concentration.  And  the  electron  density  increases 
with  temperature  in  the  freeze-out  regime  and  remains  constant  in  the  saturation  regime. 
However,  the  mobility  data  differ  in  magnitude  and  behavior,  which  further  corroborate 
the  author's  conclusion  that  the  Debye  theory  over-estimates  the  static  scieening  at  high 
donor  concentration. 


PART  V 
CONCLUDING  REMARKS 


CHAPTER  18 

SUMMARY 


Given  the  technological  importance  of  widegap  semiconductors,  it  is  highly  desired 
to  advance  the  knowledge  of  these  materials,  both  individually  and  collectively.  Despite 
the  great  diversity  in  characteristics  and  applications  of  semiconducting  materials,  there 
are  often  cases  where  systematic  understanding  is  achieved.  For  example,  the  dielectric 
bonding  theory  of  Philips  and  Van  Vechten  [6]  explains,  semi-quantitatively,  how  the 
physical  properties,  such  as  the  elastic  constants  of  CdTe,  of  tetrahedral  semiconductors 
vary  according  to  their  degree  of  ionicity.  However,  more  often  than  not,  situations 
arise  where  prior  systematic  understanding  lacks  for  the  specific  system,  such  as  for  deep 
impurities  in  phosphors.  Then  investigations  of  case  after  case  may  be  needed  to  acquire 
the  systematic  knowledge. 

This  research  is  intended,  mrough  case  studies,  to  understand  some  aspects  of  the 
widegap  semiconductors  when  the  materials  are  subjected  to  external  influences,  such 
as  quantum  confinement  and  doping.  Specifically,  modelling  of  optical  reflectivity  and 
XPS  spectrum  assists  the  implementation  of  bandstructure  calculations  and  determining 
possible  candidates  for  blue/green  electroluminescent  phosphors.  Two  salient  features  in 
an  excellent  candidate  are  required:  a  large  bandgap  (3.7  eV)  and  concentrated  density  of 
states  for  superhot  electrons  (at  3  eV).  Also,  the  modelling  of  quantum  confinement  effect 
on  the  absorption  edge  energy  of  semiconducting  crystallites  embedded  in  a  transparent 
glass  matrix  leads  to  the  realization  that  bandstructure  effect  dominates  the  contributions 
to  the.  blue  shift  of  absorption  edge  of  quantum  dots.  This  knowledge  is  useful  in  future 
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selection  of  excellent  nonlinear  optoelectronic  materials.  Furthermore,  the  investigation 
of  static  screening  effect  in  intermediately  to  heavily  doped/z-ZnSe  proves  that  the  effect  is 
indeed  important  in  determining  the  characteristics  of  the  system.  In  particular,  it  is  shown 
that  compensation  reduces  carrier  concentration  and  increases  screening  capability  of  the 
system.  Thus,  only  high-quality,  uncompensated  semiconductors  are  capable  of  efficient 
generation  of  free  carriers,  which  is  of  great  importance  in  improving  the  reliability  of 
p-n  type  injection  devices  and  implementation  of  excellent  ohmic  contacts. 

It  is  the  aim  of  basic  research  to  be  able  to  direct  future  development  by  studying 
existing  examples,  which  is  a  process  of  accumulation  of  knowledge.  Throughout  this 
research  into  several  widegap  semiconductors,  various  properties  of  them  are  computed 
from  prior  knowledge  of  their  electronic  structures  in  the  bulk  state.  This  systematic 
method  relies  upon  the  extraction  of  the  electronic  structures  from  known  experimental 
data.  This  research  fulfilled  a  tiny  loop  of  this  endless  process  of  accumulation  of  human 
knowledge,  i.e.,  application  of  the  accumulated  knowledge  to  analyze  new  data  and  in 
return  accumulation  of  new  knowledge.  This  is  what  basic  research  is  all  about.  And 
it  is  what  this  research  is  all  about. 
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APPENDIX  A 
EMPIRICAL  30  x  30  k  •  p  THEORY  [81] 

The  k  ■  p  method  perturbatively  obtains  the  electronic  structure  for  a  semiconductor 
on  the  basis  of  a  priori  knowledge  about  energy  states  at  a  high  symmetry  point  of 
k  space,  which  usually  is  the  T  point,  the  Brillouin  zone  center.  Its  basis  set  of  30 
wavef unctions  covers  an  energy  range  of  over  30  eV  and  assures  an  accurate  description 
of  the  most  relevant  valence  and  conduction  bands  in  applications.  However,  lack  of 
consideration  of  the  periodic  symmetry  in  the  reduced  Hamiltonian  lowers  its  accuracy 
away  from  the  T  point. 

Reduced  Schrodinger  Equation 

An  electron  follows  the  Schrodinger  equation  (60),  in  the  independent  electron 
approximation,  as  it  moves  in  the  periodic  lattice  of  a  semiconductor.  As  a  consequence, 
the  Bloch  theorem  allows  us  to  rewrite  Eq.  (60)  in  terms  of  the  periodic  part  of  the 
wavefunction,  un£{r).  as  equation  (65). 

H^ni{f)  =  En(k)l>n%{?)  (60) 

H=  H0  +  Vasym(f)  +  Hso  (61) 

H0=  -^-  +  V3vmm(rf)  (62) 

2  m 

ffw  =  (|)2{VV(f)xp}.c?  (63) 
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where  the.  lattice  potential,  V(r),  has  been  separated  into  a  symmetric  part,  Vsymm(f), 
and  an  antisymmetric  part,  Va3ym(f),  with  respect  to  the  interchange  of  the  anion  atom 
and  the  cation  atom  in  the  semiconductor,  and  Hso  is  called  the  spin-orbital  coupling 
term,  which  is  negligible  for  compounds  that  consist  of  elements  of  small  atomic  weight. 

*«y«  +  £ '  M*5  =  (E°  (*) "  9  v««  (65) 

Basis  Functions 

To  solve  Eq.  (65),  an  empirical  £  •  p  method  has  been  developed  [81].  This  method 
is  based  upon  group  theory,  which  allows  the  parameterization  of  the  Hamiltonian  if  a 
basis  set  of  functions  lvn  j(r),  n  =  1,  •  •  • ,  15  },  which  are  eigenfunctions  of  Ho,  is  used 
for  expanding  u  j-(r). 

HovntS(r)  =  en(o)vnz(r),      nsl,. -.,15  (66) 

where  e„(o),  n  =  1,  -■  • ,  15  are  the  energy  levels  at  the  I  point  and  the  number  of 
15  states  is  purely  the  result  of  the  group  theory  of  0/,  point  symmetry  and  a  way  of 
empirical  thinking.  It  is  enough  to  consider  only  contributions  from  the  following  states 
where  the  superscripts  /  and  u  in  state  labels  stand  for  lower  and  upper  states.  Therefore, 
only  eight  energy  parameters  are  needed  for  semiconductors  of  a  zincblende  or  diamond 
crystalline  structure. 

When  the  degrees  of  freedom  of  spins  are  considered,  the  number  of  basis  functions 
doubles,  which  is  why  the  theory  counts  30  states.  Also,  it  is  worth  mentioning  that 
the  empty  lattice  model  supplies  a  good  starting  point  to  understand  the  origin  of  these 
basis  states. 
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Table  10:  Empirical  30  x  30  k  ■  p  Theory  Basis  States 


state 

rj 

r? 

2' 

■pu 

x2' 

r! 

1  25' 

■pu 
i25' 

r« 

rw 

degeneracy 

1 

i 

1 

1 

3 

3 

3 

2 

Matrix  Elements 


Thanks  to  the  Wigner-Eckart  theorem,  the  number  of  parameters  in  the  k  ■  p  method 
becomes  manageable.  The  theorem  can  be  expressed  as 


(v['f)\Pia)\vf))  =  c\1al3;kij)P 


(67) 


where  t>[  is  the  k-th  function  that  belongs  to  the  7  representation  of  the  group  of 
the  system,  c*(7a/3;  kij)  are  the  coupling  coefficients  that  are  entirely  related  by  the 
symmetry,  and  P  depend  upon  the  physical  properties  of  the  particular  operator  P(a\ 
which  transforms  according  to  the  a-th  representation  of  the  group  of  the  system,  but 
completely  independent  of  the  quantum  numbers  i,  j,  k.  These  coupling  coefficients, 
often  called  the  Clebsch-Gordon  or  Wigner  coefficients,  can  be  found  in  reference  100 
for  the  thirty-two  point  groups. 

Applications  of  the  Wigner-Eckart  theorem  to  the  linear  momentum  operator,  the 
antisymmetric  potential,  and  the  spin-orbital  coupling  term  result  in  the  full  parameteri- 
zation of  the  Hamiltonian  H,  and  Tables  11  -13  list  the  notation  of  parameters  for  all  the 
matrix  elements.  To  express  all  the  operators  in  Eq.  (65)  in  terms  of  these  parameters, 
following  knowledge  is  necessary.  In  Oh  group,  linear  momentum  transforms  like  ^5-; 
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the  antisymmetric  potential  transforms  like  Ty',  and  the  angular  momentum  transforms 
like  Tij.  Also,  the  contribution  from  the  antisymmetric  potential  in  the  spin-orbital  cou- 
pling term  is  neglected  because  the  linear  momentum  in  the  lattice  of  a  solid  is  much  less 
than  in  atoms  and  calculations  [75]  have  justified  this.  Furthermore,  it  is  customary  to 
take  Vf  to  zero  since  the  involved  states  are  far  apart  in  energy.  Therefore,  the  number 
of  adjustable  parameters  in  the  empirical  30  x  30  k    p  theory  is  27. 

Ali  parameters  in  the  30  x  30  k  ■  p  theory  for  ZnS,  ZnSe,  CdTe  are  fitted  by  Sanders 
and  the  author  and  are  listed  in  Table  3  in  the  present  thesis. 
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APPENDIX  B 
GILAT-RAUBENHEIMER  SCHEME 


Numerical  calculations  involving  an  integration  over  the  whole  Brillouin  Zone  (BZ) 
may  be  performed  efficiently  in  two  ways:  (1)  the  special  Appoint  scheme  [101]  or 
(2)  the  Giiat-Raubenheimer  (GR)  scheme  [102,  103].  The  latter  scheme  has  been 
applied  in  the  computations  of  frequency-  or  energy-distribution  functions  of  phonons 
and  electrons,  which  ranges  from  simple  density  of  states  functions  to  generalized 
magnetic  susceptibilities,  dielectric  functions,  photoemission  energy  distributions,  etc. 
Both  schemes  allow  fine  details  of  the  distribution  functions,  such  as  the  critical  points, 
to  be  determined  as  desired.  The  present  appendix  considers  the  cubic  systems,  though 
it  can  be  generalized  to  more  complicated  systems. 

The  cubic  symmetry  reduces  the  integration  over  the  BZ  into  an  integral  involving 
only  one-forty-eighth  of  the  zone,  i.e.,  the  irreducible  Brillouin  Zone.  The  GR  scheme 
systematically  breaks  tne  irreducible  BZ,  as  depicted  in  figure  42,  into  tetrahedra  so  that, 
within  each  individual  tetrahedron,  an  analytical  interpolation  can  be  applied  to  simplify 
the  dispersion  curves.  As  a  result,  fast  numerical  implementation,  thus  the  description  of 
fine  details  of  the  distribution  functions,  becomes  possible. 

First  of  all,  the  irreducible  BZ,  labeled  as  FKWXUL  by  using  its  vertex  symbols, 
is  divided  into  three  main  tetrahedra  (itetr  =  1,2,3),  labeled  as  TLKW,  TLWIJ,  and 
FUWX.  Then,  each  of  these  main  tetrahedra  is  partitioned,  according  to  [102],  into 
smaller  tetrahedra  (t  =  1,  ■  •  •  ,n)  where  n  is  the  number  of  small  tetrahedra  in  each  of  the 
main  tetrahedra,  TLKW,  TLWU,  and  TUWX.  Within  each  of  the  subsequent  tetrahedra, 
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dispersion  curves  are  linearly  interpolated  as 


en(^J  —  tn\kicoTn)  ~H  *  f^ni^J'i  '  [^     '  ^icorn)  i 


(69) 


where  ivtx  =  l,---,4  numbers  the  four  vertices  of  the  tetrahedron.  Consequently, 
equations  (70)  and  (71)  are  derived  for  each  basic  tetrahedron  so  that  a  BZ  integral  of 
any  slowly  varying  functions  at  the  lengthscale  of  the  lattice  constant  is  converted  into 
a  summation  over  the  number  of  basic  tetrahedra. 


P?<-< 


6   E-en[k 


tetr 

( o, 


(E-e<)2 


[E-€3) 


(«3-«4)(ea-«4)(«i-«4)  ' 

(E-u)2 

{(3-(i){(2~u)((i-u)       (t3-et){c2-€i)(ei-et) 


(fi-e4)(fi- 


«3)(£i-f2)  ' 


E  <  e4 

«4  <  E  <  e3 

,       ^3  <  E  <  Co 
€2  <  E  <  6\ 

E>ei 


(70) 


where  ewtx,ivtx  =  1,  ■  ■  •  ,4  are  the  values  of  en(k)  at  the  corresponding  vertices,  v  is 
the  volume  of  the  tetrahedron,  and  £4  <  £3  <  £2  <  £1  has  been  assumed  as  well. 


(71) 


tetr 


4  / 


where  VlvU  =  tlvtx  -  E,  ivtx  -  l,"-,4  and  Dlvtx  -      J]     (Yivtx  -  Vjvtx)  with  the 

jvtx=\ 

prime  sign  denoting  that  jvtx  ^  ivtx  in  the  product.  Care  has  to  be  taken  when  treating 


the  above  expression  for  the  special  cases  when  some  of  the  Vlvtx's  are  zero  or  equal. 


143 


Figure  42:  Brillouin  Zone  of  face-center-cubic  (FCC)  and  zincblende  lattices. 


APPENDIX  C 
GENERALIZED  DEBYE  THEORY 

Basics 

At  thermodynamic  equilibrium,  we  assume  distribution  function  fi(E)  for  charged 
entity  of  type  i  at  energy  E.  The  density  of  charged  entity  of  type  i  at  energy  E  is  given  as 

m(E)  =gi{E)fi{E)  (72) 

where  gdE)  is  the  density  of  states  at  energy  E.  When  a  perturbative  electric  potential 
V(r)  is  applied,  the  density  of  charged  entity  of  type  i  used  to  at  energy  E  changes  to 

ni(E)  =  gi(E)fi[E  +  eiV{f)]  (73) 

where  the  symbol  e,  contains  the  sign  of  the  charge.  And  we  still  use  the  old  energy  to 
label  energy  states.  Therefore,  the  resultant  total  charge  density  is  given  as 

P(r-)  =  £  jdE  ti9i{E){fi[E  +  tiV{r)\  -  fi(E)}  .  (74) 

i 

The  electric  potential  satisfies  Poisson's  equation,  which  reads, 

t 
=  -  ]T  fdE  -l-g%{E){U{E  V  e,V(f)]  -  ME)} . 


(75) 


First  order  Taylor  expansion  on  the  perturbative  electric  potential  reduces  Eq.  (75)  into 
V2V(f)  =  -jV(r)  V  JdE9i(E)^l ,  (76) 
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which  allows  us  to  define  the  generalized  Debye  length  Xp  as 


A£2  =  -"  Y,[dE9i(E)- 


dE 

1 

i        j 

where  n^  is  the  concentration  of  charged  entities  of  type  i,  which  is  given  by 


(77) 


m  =  JdEgi{E)fi(E) .  (78) 


In  deducing  the  second  part  of  Eq.  (77)  we  assume  that  the  chemical  potential  Ef 
depends  upon  the  distribution  function  through  E  -  Ef,  as  always  true  for  independent 
particle  systems  where  total  energy  is  the  sum  of  energies  contributed  by  each  individual 
particles.  However,  there  is  a  catch  in  this  deduction,  that  is,  absolute  values  are  taken 
for  the  derivatives  The  reason  is  simply  because  when  real  system  is  dealt,  as  in  Part 
IV,  energy  reference  can  be  changed,  which  may  result  in  the  negativeness  of  A^2.  As  a 
note,  we  use  continuous  notation  for  energy  states  above,  which  is  only  for  convenience 
and  can  be  easily  generalized  for  discrete  cases. 

Following  Ashcroft  and  Mermin  [104],  the  relative  dielectric  constant  K\k\  can  be 


found  using 


f-\  1    Pi    ) 

K(k)  =1-—  -rk>  (79) 

v  ;  c<rv(jfc) 

where  k  is  the  wavevector  in  reciprocal  space,  p(k)  is  the  k-th  Fourier  component  of 
the  induced  charge  density  given  in  Eq.  (74),  and  v(kj  is  the  ife-th  Fourier  component 
of  the  perturbative  electric  potential  V(r)  .  The  result,  as  well-known,  is  given  below 

,(*)  =  l  +  ~j.  (80) 
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As  easily  seen,  Eq.  (77)  corresponds  to  Eq.  (43)  for  the  specific  doped  semiconductor 
system  we  discussed.  Also,  this  theory  applies  to  both  classical  and  quantum  mechanical 
systems.  The  only  condition,  which  is  sufficient,  is  that  the  total  energy  of  the  system 
is  the  sum  of  all  individual  particle  energies. 

Applications  of  Generalized  Debye  Theory 

A  direct  application  of  this  theory  on  bare  Coulomb  potential  V(r)  =  ~;  gives  the 
result  of  dressed  Coulomb  potential 

Degenerate  Electron  Gas  in  3D 
For  this  given  system,  we  have 


1  +exp  [/3{E-Ef)\ 
and  we  assume  parabolic  energy  dispersion  for  electrons,  so 

Ef<xn*-         g(Ef)  ----—,  (83) 

where  n  is  the  electron  concentration  of  the  system,  which  is  given  by 

n=  jdEg{E)f{E).  (84) 

An  approximation  is  used 
On 


r*  jdEg{E)6(E-Ef)=g(Ef) 


dEj 
so  as  to  find  an  analytical  expression  for  A#,  which  is  given  as 


(85) 


^D  oc  n    s  (86) 

and  independent  of  temperature. 
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Classical  Electron  Gas 


In  this  case,  the  distribution  function  takes  Maxwell-Boltzmann  form 

f{E)  a  exp  [*.fi{E  -  Bf)] 

and  no  matter  what  energy  dispersion  for  electrons  is  we  have 
dn  d 


dEf  "  dEf 


IdE  g{E)Cexp  [-/?(£  -  Ef)}  =  f3n 


so  that 


(87) 


(88) 


\j)  ex  n"'2T2 


(89) 


which  is  quite  different  from  the  previous  case  in  its  dependence  on  both  the  temperature 
and  carrier  concentration. 
Singly  Charged  Donors 


In  the  present  case,  the  number  of  neutral  donors  is 

Nd_ 

1  +  le-P(E4+Ef) 


N*d 


according  to  Eq.    (41)  and  we  have 


A", 


N+  =  — 
d   "  1  -r  2eHE*+Ef)  ' 


Therefore,  it  is  easy  to  find  that 


1  8Ef  '  ""  P  Nd 


(90) 


(91-) 


(92) 


and  owing  to  the  symmetry  in  their  relationships,  it  is  trivial  to  find  the  solution  for  the 
system  of  singly  charged  acceptors,  for  which  we  list  the  result 


Ma-  Na~{Na-N-) 

dEf  l  '"  P  Na 


(93) 
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